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ABSTK/CT 


The  Application  of  reliability  analysis  methods  to  the  estimation  of 
probable  aircraft  structural  fatigue  performance  was  investigated,  use 
was  made  of  order  statistics  to  establish  the  means  of  assessing  the 
fatigue  performance  reliability  of  a  fleet  or  number  of  'atlgue-exposed 
details.  A  reliability  analysis  plan  for  application  to  aluminum  alloy 
structural  fatigue  performance  was  developed  and  compared  ««rh  the 
current  f lxed-scatter-fectoi  procedure  for  determining  the  safe  life 
of  a  structural  detail.  Both  ti.e  two-parameter  Veibull  distribution 
and  the  log-normal  distribution  with  empirically  defined  sh.tpe  parameters 
were  used  to  make  the  reliability  plan  tractable  as  compared  to  a  dis¬ 
tribution-free  approach .  Max tmum- likelihood  estimators,  including  one 
th&t  considers  only  the  first  two-ordered  failures,  were  employed  to 
examine  the  many  variables  that  might  influence  fatigue  scatter,  to 
qualify  fatigue  date  that  represented  aluminum  structural  scatter,  and 
to  establish  shape  parameter  values  that  typified  structural  fatigue 
scatter .  The  sampling  distributions  of  these  estimators  were  required 
to  work  the  problem  and  were  calculate  1  by  means  of  existing  theory  or 
Monte-Carlo  simulation.  More  than  2,000  groups  of  fatigue  performance 
data  were  collected,  analyzed,  and  used  to  demonstrate  the  feasibility 
of  establishing  a  shape-parameter  value.  Based  on  this  estimate , 
scatter  factors  have  been  generated  to  account  for  the  penalty  of 
limited  input  Information,  the  degree  of  desired  reliability,  and  the 
sise  of  the  exposed  fleet.  Using  these  factors,  the  possible  effects 
of  the  reliability  analysis  on  structural  weight,  payload,  or  range 
wers  explored  for  a  let-engined  military  tanker /transport-type  airplane. 


This  abstract  is  subject  to  special  export  controls,  and  each 
transmittal  to  foreign  government  or  foreign  nationals  may  be 
made  only  with  prior  approval  of  tbs  Air  Force  Materials 
Laboratory,  Wright -Patterson  APB,  Ohio  45433. 
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F  •  when  used  alone,  probability  of  failure  event  or, 

when  used  with  s  random  variable,  a  symbol  of  any 
time-to-fal.'ure  distribution  function 

.  U 

»  probability  of  failure  of  the  component  destined  to 

or  have  the  Mc**  lowest  fatigue  life  of  I*  Independent 

exposed  components 

■  the  1th  ordered  censored  test  observation,  where  the 
test  or  exposure  has  been  terminated  for  reasons  other 
than  failure 

-  population  average  of  log  failure  tlmaa,  when  the 
population  distribution  is  not  necessarily  assumsd 
• o  be  log-normal 
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h  ”  hazard  function;  also  known  as  the  intensity  or 

failure  rate  function  or  as  tne  force  of  mortality 

<  -  1/a,  the  We. bull  parameter  that  la  actually  eetlaatad 

in  this  report,  rather  than  a  itself 

k  -  the  standardized  normal  deviate  corresponding  to  boom 

p  probability  nuaber,  p.  More  than  one  definition  has 

been  discovered  in  the  literature;  in  this  report, 

k  is  defined  bv 
P 


M 


h‘ 


n 

n ^  or  k 


l 


exp  (-£/  /j)  df. 


P 


-  any  unknown  parameter 

■  likelihood  function 

■>  log-nonul  distribution  scale  paraaeter  or  nedian 
life 

«  order  nuaber  given  to  the  ^  lowest  fatigue  life  in 
a  fleet 

-  fleet  size;  the  nuaber  of  noainally  identical, 
Independent  structural  coaponents  exposed  to  one 
fatigue  envlronaent 

-  test  ssaple  size 

■  nuaber  of  failure  observations  in  the  test  ssaple 


Q 


n 


R 


*M;N 
or  R 


■  n-ttj ,  the  nuaber  of  observations  not  Involving 
failures 

■  coefficient,  depending  on  n,  proportional  tc  the 
variance  of  an  estimator's  aaapling  distribution 

•  1-F;  when  used  alone,  R  is  the  single,  randoaly 
chosen  fleet  aeaber’s  reliability  or  probability  of 
no  failure.  When  used  with  a  randoa  variable,  it.  is 
s  syabol  denoting  any  reliability  function. 

-  1-F,  the  reliability  of  the  structural  coa^nwnt 

destined  to  have  the  Hth  lowest  fatigue  life  in  a 
fleet  of  N  identically  exposed  independent  coaponsnts 


c  or 
a (leg  X) 


2 

a 


log-noraal  shape  paraateter 

log-noraal  paraaeter  actually  est lasted  lr  this 
report,  rather  tnan  o  itself 
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-  sample  variance  of  the  logarithm  of  failure  observe- 
tions  only  f roe  a  test  sample 

-  total  scatter  factor  (  -  S  Ss  •  Su  ) ;  defined 

nf  R  kr 

aa  the  ratio  of  the  beat  estimate  of  the  scale 
parameter  (characteristic  life  or  average  log  life) 
to  the  final  (interval)  estiute  of  safe  life 

-  factor  to  account  for  finite  siae  of  the  full-scale 
fatigue  test  sampling 

-  factor  required  to  provide  a  desired  level  of 
reliability  for  a  fleet  of  just  one  airplane , 
component ,  or  detail 

•  factor  to  account  for  fleet  siae  exposed;  double 
subscript  is  used  where  this  factor  also  depends 
on  the  preassigned  reliability 

■  S«  *  S  ,  factor  required  to  provide  reliability  R 

*  R 

to  a  fleet  of  else  one  or,  equivalently  to  provide  a 
lower  reliability  R  to  all  or  moat  of  a  fleet  of  size 

H 

-  a  discrete  random  variable  denoting  the  numbe'  of 
failures  occurring  in  at.  exposed  fleet 

-  the  statistic  a/a  (see  definitions  of  accent  marks  that 
follow  ) 

V  a 

-  the  statistic  (8/6) 

»  the  statistic  (8/B)a 

■  symbol  of  the  chi-squared  distribution 

■  the  y-fractile  of  the  chi -squared  distribution  with  f 

degrees  ct  freedom 

»  a  continuous  tiae-to-f allure  random  variable  usually 
referring  to  teste  useful  in  estimating  fatigue 
scatter  only 

*  log  X 

»  e  time-to-f allure  variable  usually  referring  to  the 
fatigue  performance  of  one  structure  or  structural 
component  in  service  or  to  fatigue  results  from  full- 
scale  teats  useful  in  estimating  mean  life 
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a  tine-co-failure  random  variable  denoting  the  Xth 
lowest  fatigue  life  In  a  fleet  service  exposure 
of  else  N 


safe  life,  l-  represents  the  (1-R)  fractile  of  the 

distribution  of  Z;  YR  represents  the  (1-R)  fractile 

of  the  distribution  of  Y.  ("Life"  is  defined  aa  the 
tine  to  initial  appearance  of  fatigue  cracking  et  a 
detail  or  in  a  structure.) 

NOTE:  The  standard  practice  of  using  uppercase  letters  to 

denote  randan  variables  and  corresponding  lowercase 
letters  to  denote  miner  leal  values  the  variables  nay 
take  is  followed,  in  nost  instances,  in  this  report. 

ACCENT  MARKS  DENOTING  ESTIMATION 

(a)  Bar ,  as  in  or  A ,  nay  denote  an  average  or  any  sort 
of  p<'lnt  aatiaator.  ’ t  is  to  be  noted  that  the  bars  in 
F  and  R  do  not  have  this  neanl "g  (see  Par.  2,  "Symbols , " 
above)  . 

(b)  Boat .  as  in  TR,  denotes  e  one-sided  bound  of  a  confi¬ 
dence  interval  estimate.  This  accent  nark  nay  be 
replaced  by  or  used  in  conjunction  with  the  subscript  y 

A 

as  in  6  or  6  . 

Y  Y 

ic)  Hat ,  as  in  o,  symbolises  the  maxinus-llkslihood  point 

estimate 

(d)  Double  hut,  as  in  a  represents  an  estimator  (usually  the 
HUE)  that  is  free  or  has  been  freed  of  bias. 

(e)  Prime,  as  in  denotes  an  estimator  taken  from  acre  than 
one  sample. 

(f)  Sjulmlt,  as  in  f,  stands  for  an  unbiased  two- ordered- 
statistic  estimator  of  a  shape  parameter. 

v 

(g)  Vaa ,  as  in  8,  denotes  a  scale-parameter  MLR  that  aust 
depend  on  an  estimated  value  rather  than  on  a  true 
shape-par ase ter  value. 


SECTION  I 


INTRODUCTION 


The  structural  reliability  of  an  airplane  ta  dapanoant  on  tha 
likelihood  of  a  fatigue  failure  or  a  static  failure,  or  a  combination 
of  both.  At  the  present  stage  In  the  evaluation  of  the  aluminum  alloy 
airframe  the  expectancy  of  fatigue  crack  initiation  can  be  of  a  higher 
order  than  that  of  a  static  failure.  On  the  other  hand,  it  should  be 
realised  that  with  a  well-designed,  fail-safe  structure,  subjected  to 
regular  and  correctly  scheduled  inspections,  the  fatigue  damage  initia¬ 
tion  incident  should  not  progress  beyond  the  detectable  crack  atage  to 
the  level  tdiere  a  structural  failure  la  imminent.  The  aircraft  designer 
still  has  a  real  requirement  to  assess  the  expected  time  at  which  cracks 
may  Initiate  on  a  particular  structure.  In  order  to  Insure  a  sufficiency 
of  trouble-free  perfcmunce .  this  report  la  primarily  concerned  with 
the  development  of  a  reliability  analytic  method  for  application  to  thia 
fatigue  problem  as  affactad  by  safsty  or  maintenance  requirements. 

Some  clarification  of  the  term  "reliability  analysis"  la  necessary 
as  its  application  to  atructuraa  la  not  tha  usual  one  of  camponant 
reliability  leading  directly  to  the  system's  reliability  or  to  total 
rsliability  through  an  accounting  of  tha  interdependence  of  components, 
fault-tree  analyses,  or  ocher  such  procedures.  Structural  reliability 
as  envisaged  here  la  concerned  with  a  specific  structural  detail  and  is 
heavily  dependent  on  the  tasting  of  full-acale  structural  details  or 
prlauary  structural  components .  For  axampla,  Ralthby  (1)  has  shown  that 
thia  form  of  tasting  is  easantlal  to  obtain  an  estimate  of  available 
fatigue  performance  that  bears  any  res ambiance  to  that  actually  achieved. 
It  has  also  bean  observed  that  structural  fatigue  problems  are  fre¬ 
quently  initiated  by  fretting,  a  process  that,  is  not  reproduced  in  the 
%re  11 -designed  simple  laboratory  specimen  and  its  teat  program.  Like¬ 
wise,  estimates  of  the  mean  performance  baaed  on  such  simple  tests 
would  probably  be  larger  than  that  demonstrated  by  a  structure  in 
service.  The  significant  goal  of  a  reliability  analysis  is  to  determine 
the  probable  time  to  the  Initial  appear ancea  of  fatigue  daamge  in 
structures  or  details,  a  task  even  more  obviously  beyond  the  scope  of 
the  simple  specimen. 

Other  ’’sal  hasards — bayoud  those  of  the  teat  specimen — confront  the 
designer.  Load  environment,  physical  environment,  and  actual  local 
stress  response  of  a  detail  to  its  loading  environmant  all  lend  a  con¬ 
founding  element  to  the  task.  Gross  misrepresentation  of  airplane 
fatigue  performance  through  inadequate  fatigue  test  design  and  interpre¬ 
tation  will  generally  lead  to  poor  rallabillty  estimates.  However,  It 
le  as  meaningless  to  aaeociata  these  errors  with  the  reliability  analysis 
as  It  would  be  to  regard  the  errors  in  the  assumed  functional  form  of 
the  fatigue  life  distribution  as  weaknesses  in  Che  structural  analyses. 
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Currant  analytical  practice  (Ref .  2)  for  allltary  aircraft  for 
insuring  an  adequate  or  aafe  level  of  fatigue  performance  by  aircraft 
■tructurea  doe*  contain  acme  of  the  characterletice  of  a  reliability 
plan.  Namely,  a  structural  detail,  component,  or  complete  airframe  la 
fatigue- tea ted  to  a  defined  and  representative  loading  environment. 
Additionally,  a  goal  fa  selected  at  some  multiple  of  the  defined  and 
expected  eervlce  life  of  the  structure.  The  applied  eafety  factor  la 
usually  aalactad  quit#  arbitrarily  aithar  on  the  baaia  of  judgment 
or  with  guidance  from  soma  possible  consideration  of  confidanca  and 
probability  levels  of  statistical  analysts. 

In  past  assessments  of  fatigue  scatter,  the  log-normal  distribution 
has  oftan  been  used .  Attsntion  has  also  been  given  to  the  Welbull-type 
distribution  for  defining  fatigue  variability  and  predicting  potential 
fatigue  performance.  The  important  consideration  is  that  a  level  of 
fatigue  performance,  statistically  defined  or  otherwise  considered 
acceptable,  is  set  for  the  structural  detail.  A  vital  part  in  this 
assessment  has  been  the  use  of  mean  or  median  fatigue  performance  for 
the  detail.  However,  these  part  analyses  have  not  considered  the  likely 
appearance  of  early  or  first  fatlgus  damage  during  the  exposure  of  a 
quantity  of  details  to  the  fatigue  environment.  The  expected  time  to 
flret  failure,  rather  than  either  a  probable  level  or  the  average  per¬ 
formance  with  Its  implied  high  number  of  failures  (about  half  of  the 
fleet,  prior  to  the  design  goal),  has  been  proposed  by  Freudenthal  (3) 
as  an  important  design  paremeter  for  structural  reliability.  The  early 
failure  is  also  a  serious  economic  consideration  that  could  possibly 
influence  future  decisions;  consequently,  this  added  factor  must  be 
investigated  and  carefully  weighed  by  the  design  engineer.  It  is  known 
that  if  fstigue  life  comes  from  a  population  whose  probability  density 
is  of  a  specified  limited  type,  Che  extreme  values  (initial  failures) 
of  an  indefinitely  large  sample  will  have  a  Welbull  distribution.  For 
the  problem  of  the  clrplane  fleet,  where  a  substantial  number  of 
structural  details  are  involved,  the  use  of  this  asymptotic  form  oet>« 
appropriate.  Some  recent  papers  describing  or  using  this  approach  are 
available  In  tha  literature  (A  through  9),  and  this  report  is  r- .  exten¬ 
sion  or  developawnt  of  this  work,  with  special  emphasis  on  th~.  estima¬ 
tion  procedure  part  of  the  problem. 

Section  II  presente  a  complete  outline  of  the  proposed  reliability 
analysis  and  makes  comparisons  with  the  existing  procedure  for  determin¬ 
ing  the  safe  life  of  an  aircraft  structure.  The  plan  has  been  divided 
into  five  steps.  The  first  four  steps  describe  the  assumptions, 
estimating  procedure,  and  the  influence  of  sample  size.  The  final  step 
defines  the  influence  of  fleet  or  number  of  detail  exposed  to  the 
fatigue  environment. 

Section  III  contains  discussions  on  the  processing  techniques  used 
on  the  mesa  of  collected  date,  the  procedutee  developed  for  coping  with 
problems  that  arose  during  processing ,  end  the  results  that  were  obtained 
from  the  evaluated  fatigue  performance  date.  t  dissertation  is  presented 
on  a  revised  scatter-factor  approach  for  ca)  ■■  elating  the  safe  life,  a a 
defined  in  Sec.  II,  of  an  aluminum  structure,  the  possible  performance 
penalties  suffered  by  adoption  of  this  proceed  method  are  also  discussed . 
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Section  IV  presents  «n  outline  of  the  Analytical  groundwork 
needed  for  estimating  the  distribution  shape  parasMters  from  the  ease 
of  collected  fatigue  data.  The  techniques  for  obtaining  blaa  and 
variance  values,  the  net hod  of  weighting  estimates  frost  groups  of 
varying  sample  slse,  and  the  procedure  for  obtaining  interval  esti¬ 
mates  are  presented. 

Section  V  sum arises  all  the  additional  fonsulatlona  that  were 
used  for  estiaating  s  safe  life  for  an  aircraft  structure,  which  include 
the  techniques  for  calculating  the  scale  parameters  and  the  ao-called 
scatter  factors. 

Section  VI  presents  a  list  of  the  conclusions  aade,  together  with 
s  few  suggested  recommendations , 

Appendix  I  contains  proofs  for  the  important  theorems  used  in  the 
analysis,  end  aoaie  aspects  of  the  Monte-Carlo  techniques  that  were 
adopted  are  discussed. 

Appendix  II  lists  all  the  source  references  and  tabulates  the 
salient  features ,  including  the  point  estlauates  of  the  shape  and  scale 
parameters,  of  the  collected  fatigue  performance  data. 

Finally,  Appendix  III  lists  the  cyclic  life  for  every  individual 
specimen  assessed  in  this  report . 
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SECTION  II 


OUTLINE  OF  THE  RELIABILITY  ANALYSIS 


The  basic  objective  of  this  investigation  is  the  application  of 
reliability  analysis  to  aircraft  structural  fatigue  performance,  in 
terms  of  the  first-  and  second-ordered  times-to-fat igue-daaage  initia¬ 
tion.  Statistical  definition  of  fatigue  variability  is  essential  to 
the  achievement  of  this  objective. 

The  generality,  ucfinttion,  and  input  data  requirements  of 
distribution-free  Methods  may  make  their  application  impractical  and 
economically  forbidding.  Furthermore,  rationale  juat  seems  to  point 
toward  a  more  tractable  definition  of  fatigue  variability  by  some 
reasonable  and  determinable  distribution.  The  two-parameter  log-normal 
distribution,  bv  virtue  of  its  easy  transf ormation  from  the  well-known 
normal  or  Gauaalan  distribution,  has  been  a  useful  engineering  tool  for 
evaluating  reliability  in  fatigue  performance.  The  two-parameter 
Welbull  distribution  function  has  also  found  an  application  to  fatigue 
variability  studies.  Presuming  knowledge  of  the  shape  parameters  of 
thase  distributions ,  statistical  analysis  procedures  that  account  for 
sample  size  are  readily  developed.  Using  these  procedures,  the  analyst 
can  calculate  interval  estimates  of  the  true  fatigue  life  distribution 
at  any  desired  level  of  confidence.  The  additional  tolerance  for 
probability  of  failure  of  any  one  detail  la  then  also  definable.  Order 
statistics  carry  the  fatigue  reliability  analysis  a  step  forward  by 
focusing  on  the  likelihood  of  the  first  or  the  second  failure  in  an 
exposed  fleet  or  in  groups  of  details.  In  any  cate,  development  of  a 
reliability  analysis  as  either  a  probable  levee  of  fatigue  performance 
or,  even  more  detailed,  in  terms  of  an  ordered  statistic  has  consider¬ 
able  complication  over  the  simple  exercise  of  applying  a  scatter  factor 
to  an  estimated  average  or  median  fatigue  performance. 

7he  main  effort  of  this  research  project  has  been  devoted  to  the 
development  of  tha  reliability  analysis  tools  that  can  be  applied  to 
the  aircraft  structural  fatigue  performance  ta^k.  Attempts  have  been 
made  to  fully  utilize  all  applicable  fatigue  data  and  to  employ  only 
mathematically  precise  estimation  procedures  in  the  analysis  and 
dstaction  of  important  variables  and  vagaries  of  structural  fatigue 
performance. 

The  analysis  does  not  account  for  all  unknown  or  inadvertent 
systematic  errors  that  may  cause  the  full-scale  fatigue  test  result 
and/or  the  ensuing  life  calculations  to  grossly  misrepresent  the  actual 
aircraft  fatigue  performance.  Such  errors  could  arise  from  the  assumed 
load  environment ,  customer  usage,  and  cumulative  damage  theory  aa  veil 
as  from  the  specimen  design  and  fatigue  teat  setup.  In. all  that  follows, 
it  is  assumed  that  the  life  estimate  arising  from  full-scale  teats  ia 
drawn  from  the  same  statistical  population  as  are  Che  corresponding 
aircraft's  structural  lives  and  that  this  population  distribution 
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represents  fatigue  strength  variability  only.  In  analytic  strength 
variability,  specific  attention  hae  been  given  to  accounting  for  the 
sample  six*  of  full-scale  developmental  fatigue  tests  and  for  the 
eventual  alee  of  the  fleet  to  be  exposed  m  s  defined  fatigue  environ¬ 
ment.  This  ha#  been  done  through  the  application  of  statistical  Infer¬ 
ence  theory  and  of  order  statistics  concepts  respectively. 

The  need  for  this  accounting  Is  clear  when  the  current  approach  to 
create  reliability,  represented  by  Fig.  1,  le  compared  to  the  proposed, 
more  comprehensive  approach  outlined  and  described  schematically  in 
Fig.  2. 

In  the  current  approach,  the  symbol  Y  in  Fig.  1  is  a  measure  of 
the  uvei age  performance  of  developmental  fatigue  , c o fc  '"cc Isctu  and  la 
usually  taken  to  be  the  mean  of  log  t tmeu  or  cvcles  fo  failure.  This 
average  performance  is  then  divided  by  a  so-called  scatter  factor, 
whoa#  value  may  range  from  two  to  four,  to  obtain  Ys.  which  Is  tradition¬ 
ally  known  as  the  "safe  life,"  a  life  at  which  the  probability  of  the 
Initial  appearance  of  fatigue  cracking  at  a  detail  is  sufficiently  low. 
These  values  of  the  scatter  factor  have  been  subject  to  Interpretation. 
Hopefully,  these  factors  would  provide  an  estimate  of  the  service  life 
prior  to  initiation  of  fatigue  c racking.  In  other  circumstances,  thev 
are  credited  with  providing  a  fatigue  performance  with  reasonable  main¬ 
tenance  throughout  the  economical  life  of  the  structure.  Fall-safe 
design  concepts--wlch  their  damage  containment  capability — have  con¬ 
sidered  that  the  magnitude  of  the  factor  can  be  reduced  over  that 
required  for  non-fail-safe  structure.  Seemingly,  the  scatter-factor 
value  chosen  is  independent  of  the  number  of  details  to  be  exposed  in 
service  and  of  the  sample  sire  of  the  full-scale  airplane  test(»). 

In  contrast,  the  difference  between  mean  test  performance  and  sub¬ 
sequently  calculated  probable  or  safe  life  is,  in  the  proposed  approach, 
strongly  dependent  on  sample  and  fleet  sire.  Recognition  of  the  rela¬ 
tive  Importance  of  the  structural  detail,  including  the  degree  of 
redundancy,  may  still  be  accomplished  using  the  plan  outlined  in  Fig.  2. 
However,  in  this  new  approach,  an  acceptable  level  of  reliability  must 
be  specified.  Scatter  factor  may  then  be  calculated;  to  specify  scatter 
factor  a  priori  is  equivalent  to  specifying  a  reliability  that  is  almost 
certainly  not  of  the  desired  level, 

As  the  first  of  five  major  steps  in  applying  statistical  techniques 
it  wss  decided  to  consider  the  untruncated  two-paraaetcr  Weibull  distri¬ 
bution  aa  a  random  t iae-to-f allure  model  because  of  certain  desirable 
mathematical  properties  It  possesses  (these  are  discussed  In  Sec.  III). 
That  la,  the  fatigue  life  Y  is  assumed  to  be  a  random  variable  with 
distribution  given  by 

F(y)  *  1-exp  {-(y/8)a)  for  unknown  a, 8  >  0  < 1 1  -  1 ) 

The  traditionally  cited  log-normal  model  was  also  applied  to  the 
problem  in  a  less  complete  manner  due  to  difficulties  In  working  with 
estimates  of  its  two  parameters  from  incomplete  or  censored  samples. 
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These  ere  samples  In  which  Che  testing  of  one  or  more  specimens  is 
regarded  es  being  terminated  for  reasons  other  than  failure. 

At  this  time,  no  attempt  has  been  made  to  test  the  validity  of  any 
erne  distribution  or  even  to  reject  one  model  for  another.  However, 
somewhat  different  conclusions  may  be  reached  from  applying  different 
distribution  functions  to  the  same  data.  Of  course,  the  necessity  to 
generate  very  large  samples  of  controlled  fatigue  data  may  preclude 
adequate  resolution  of  this  problem.  To  emphasize  the  need  for  making 
realistic  assumptions  about  the  distribution  of  fatigue  lives,  the 
Tchebycheff  inequality  will  also  be  applied  to  the  problem  to  indicate 
a  possible  conclusion  when  no  assumption  is  made  about  the  functional 
form  of  the  true  time-to-failure  distribution. 

Given  reasonable  accuracy  in  the  chosen  model,  the  second  step  is 
the  most  important,  necessary,  and  critical  step  in  the  analysis.  It 
was  assumed  that,  for  a  wide  variety  of  structure  and  structural  simu¬ 
lation  configurations  and  types  of  loading,  the  variability  or  disper¬ 
sion  of  the  logarithms  of  times-to-failure  is  nearly  constant.  More 
specifically,  with  reference  to  the  distribution  function  itself,  the 
assumption  is  that  the  scatter-controlling  or  shape  parameter  of  the 
assumed  Welbull  (or  log-normal)  model  has  a  unique  value  for  this  same 
9et  of  conditions.  This  unknown  value  of  the  Weibull  shape  a  has 
been  estimated  from  a  collection  of  more  than  a  thousand  applicable 
samples  of  two  or  more  fatigue  specimens  each,  and  the  errors  in  the 
uniqueness  assumptions  are  considered  to  be  small  enough  for  practical 
application.  A  large  proportion  of  Che  later  sections  of  this  report 
in  devoted  to  discussing  and  justifying  this  important  conclusion. 

With  the  shape  parameter  a  accurately  estimated  and  subsequently 
regarded  as  being  fixed,  we  need  only  obtain  from  full-scale  testing  an 
accurate  estimate,  F,  of  the  scale  parameter  B  to  complete  the  third 
stage  of  the  plan.  This  step  establishes  Curve  A  In  Fig.  2  and  pro¬ 
duces  the  best  point  estimate  possible  of  the  complete  distribution  of 
fatigue  lives.  It  is  from  this  curve,  with  equation 

P(y)  -  1-exo  [-(y/B)aJ  (II-2) 

that  the  most  likely  probability  of  failure  at  a  given  fatigue  life  is 
es  tlmated . 

It  is  standard  procedure  in  any  engineering  discipline  to  specify 
allowables  as  values  somewhat  less  than  the  most  likely  estimates.  One 
sensible  way  to  do  this  was  chosen  as  step  four,  where  a  confidence 
interval  estimate  #  of  B  is  made  that  accounts  for  the  finite  size 
of  the  full-scale  fatigue  sampling  so  that 

P(^y>  B)  -  y  (II-3) 
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The  confidence  coefficient,  y  i#  generally  specified  et  ease 
acceptable  level.  The  estimate  ty  along  with  the  "known"  shape 
p*remeter  a  completely  determines  Curve  B,  a  plot  of  the  interval 
estimate  of  the  complete  distribution,  with  aquation 


f  (y)  -  l~e*p  [-(y/fl  )a] 

y  ‘  Y 

It  follows  immediately  from  Eqe.  ( 1 1  - 3 )  and  (II-4)  that 


( II-4) 


P  f^(y)  *  F(y )  1  -  y  (11-5) 

In  other  words,  with  probability  y,  the  estimated  failure  probability 
(reliability)  from  Eq .  (II-4)  will  be  greater  than  (lees  than)  the 
unknown  true  failure  probability  (reliability). 

As  s  structural  design  specification,  the  type  of  confidence  state¬ 
ment  given  above  Is  certainly  not  without  precedent.  For  example ,  the 
"A"  and  "B"  values  of  ultimate  strength  of  structural  alloys  quoted  in 
MIL-HDBK-SA  (Ref.  10)  are  baaed  on  a  confidence  interval  estimate  analy¬ 
sis  that  differs  from  the  sf  or  Mentioned  one  in  only  two  respects.  The 
first  difference  is  that  the  random  variable  ultimate  static  strength 
was  assumed  to  have  r  two-parameter  normal  distribution.  Ths  second 
difference,  one  of  greater  practical  importance,  is  that  both  param¬ 
eters  were  assumed  to  be  unknown  in  Ref.  10  and  ware  estimated  from  one 
large  sampling,  of  at  least  100  observations,  of  tensile  specimens.  As 
outlined  before,  the  shape  parameter  vaa,  in  this  study,  asetaMd  to  ba 
coamon  to  the  actual  aircraft  aa  wall  as  to  the  large  collection  of  data 
from  which  it  waa  estimated.  Some  assumption  regarding  prior  informa¬ 
tion  about  th._  degree  of  fatigue  scatter,  as  measured  bv  the  shape- 
parameter  estimates,  is  obviously  required  in  this  study  becausa  the 
full-scale  fatigue  teat  aaaple  slat  typically  does  net  exceed  one  or  two 
observations .  Two  parameters  cannot,  of  course,  ba  estimated  at  all 
from  a  sample  of  one  and  cannot  be  estimated  accurately  from  a  sample 
of  two  observations . 


The  first  four  stages  of  the  reliability  plan  have  lad  to  Curve  B, 
which  is  simplv  a  statement,  at  some  confidence  level,  of  the  reliability 
of  an  aircraft  or  aircraft  structural  detail  (depending  on  the  nature 
of  the  full-scale  teat),  chosen  randomly  from  s  universe  of  nominally 
identical  counterparts ,  as  a  function  of  time  of  exposure  to  the  load¬ 
ing  environment.  If  the  reliability  R  of  an  arbitrary  fleet  member  is 
judged  to  be  a  aultabla  maaaure  of  safety  or  of  freedom  from  repair,  the 
plan  terminates  upon  administrative  specification  of  an  acceptable  value 
of  R  and  upon  reference  to  Curve  B  to  certify  a  corresponding  "safe 
life." 


In  many  or  possibly  in  most  instances,  the  arbitrary  specification 
of  an  acceptable  value  of  R  Is  not  a  sufficient  basis  for  certifying 
safe  life  for  a  number  of  details  or  i  fleet  of  airplanes.  With  regard 
tc  the  safe  or  repair-free  operation  of  a  flaet,  both  customar  and 


vendor  usually  expect  that  the  entire  or  neatly  the  entire  fleet 
demonstrate  acceptable  reliability.  Precise  definition  le  one  step 
toward  this  goal. 

With  reference  to  the  proverbial  chain  that  la  no  stronger  than  ite 
weakest  link,  it  la  only  a  natter  of  definition  to  equate  the  relia¬ 
bility  of  the  majority  of  the  fleet  with  the  reliability  of  one  of  the 
weakest  of  ite  (assumed)  Independent  nominally  identical  members.  That 
ie,  the  reliability  of  the  entire  fleet  le  not  different  from  the  relia¬ 
bility  of  its  weakest  member.  In  general,  the  reliability  of  the  vest 
majority  of  the  fleet,  considered  as  a  package  and  excluding  only  the 
weakeat  H- i  members .  is  defined  as  the  reliability  of  the  M1*1  weakest 
fleet  member. 


The  random  variable  ZM;N,  used  to  denote  the  fatigue  life  of  the 
detail  destined  to  have  the  Nth  lowest  of  N  fatigue  lives,  la  an  example 
of  an  order  etetlatic.  It  has  its  own  statistical  distribution.  This 
distribution  is  completely  determined  by  the  randomly  chosen  member 'a 
fatigue  life  distribution.  Note  chat  the  two  distributions  are  never 
equal  except  in  the  trivial  case  of  M  ■  H  ■  1. 

Fortunately,  the  precise  calculation  of  the  order  atetistic'a 
distribution  from  any  assumed,  estimated,  or  known  parent  distribution 
function  is  a  simple  matter  whenever  the  trials  are  independent,  how¬ 
ever  complex  the  parent  function.  The  derivation  of  the  cumulative 
binomial  distribution,  which  msy  always  b#  used  to  asp  ths  parent  into 
the  order  statistic ‘a  distribution,  is  given  in  many  standard  statis¬ 
tical  texts  (e.g.  kef.  11).  Xa  Fig.  3  this  summed  binomial  formula  le 
plotted  for  several  combinations  of  M  and  N.  Tha  raeultant  curves 
are  worth  eziwinatlon  by  ell  readers,  but  they  will  be  especially  en¬ 
lightening  to  those  not  fanllis.  with  the  behavior  of  order  statistics, 
since  they  point  out  significant  trends  that  ere  perhaps  not  obvious. 

The  most  striking  example#  erlsa  from  the  way  the  curve  slopes  Increase 
with  M,  leading  to  such  surprising  facta  as: 

1)  Greater  fleet  member  reliability  R  and  tnerefore  larger 
scatter  factors  and  more  weight  on  the  airplane  results  from 

a  specified  reliability  of  99Z  for  the  weakeat  of  10  than  from 
the  a (  *  specification  for  the  second  weakest  of  100.  The 
revers..  ie  true  for  leaser  reliabilities  such  as  f  »  0.90. 

2)  In  a  similar  vein,  if  the  reliability  of  the  single,  randomly 
chosen  fleet  member  la  greater  then  0.97,  it  is  more  likely  to 
experience  one  or  more  failures  in  a  fleet  of  size  4  than  it 
ia  to  experience  two  or  more  failures  in  a  fleet  of  20.  Por 

R  <  0.97,  this  trend  la  reversed. 

This  trend,  which  depends  vitally  on  the  independence  of  fleet  members 
but  which  in  nc  way  depends  on  the  assumed  failure  distribution  could 
be  Interpreted  by  the  engineer  or  administrator  in  a  practical  manner. 

It  points  out  the  fact  that,  except  for  cracks  or  failures  that  are 
"intolerable,**  It  would  be  quite  worthwhile  to  be  able  to  specify  relia¬ 
bility  for  the  second  or  third  weakest  fleet  seaber  rather  than  to 
specify  a  similar  reliability  level  for  the  weakest  fleet  member. 
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Returning  to  the  final  itagt  of  the  reliability  plan  as  outlined 
in  Fig.  2,  it  la  proposed  that  an  acceptable  reliability  Itwl  %!}* 
b«  specified  f row  which  the  cueulatlv*  blnonlal  distribution  can  be  need 
to  calculate  the  correa ponding  ‘general"  reliability  t.  The  certifiable 
safe  life  nay  then  be  found  fron  Curve  B  to  coaq»lata  the  eatarcla*. 

In  sunnary,  the  object  of  the  overall  atatiatlcal  aanlyaia  la  to 
obtain  an  accurate  and  reliable  eat  legate  Iff  of  the  "true”  design  lift 
Zf ,  where  Zf  eatlaflea  the  probabilistic  atateaent: 

The  proportion  of  the  failure  tinea  Z  that  will  exceed  life 

Zf  ia  equal  to  whcte  Z  denote#  the  Mth  failure  tine  In 

each  flaet  of  a  univerae  of  flaata  of  a  lea  N. 

In  equation  f one  thia  le  expressed  ae 

P(Z  *  2j)  -  l  (II-*) 

Since  Zf  can  only  ba  aatlnated  froa  eanpla  data  wa  can,  at  bant, 
obtain  an  estimate  ¥f,  which  ia  cooaervatlva  with  confidence  y,  euch 
that 


P(^  <  If )  -  y  (II-7) 

Specification  of  levala  of  reliability  and  confidence,  the  fatigue 
taat  eanpla  else,  and  the  fleet  exposure  nust  r ana in  a  Joint  Managerial 
effort  between  cue toner  and  vandor.  It  ia  recognised  that,  aa  a  prac¬ 
tical  natter,  thaee  epecif i cation*  will  depend  pertly  on  their  resulting 
calculated  ecatter  factors  and  ultlaataly  on  structural  weight.  The 
relationship  of  structural  weight  and  level  of  reliability  ia  explored 
in  Sec.  III. 
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SECTION  III 


DISCUSSION 


1.  DATA  EVALUATION  PROCEDURES 
a.  General  Procedure 

Tha  choice  of  the  Weibull  etat let  leal  nodal  can  ba  justified  on 
the  baeie  that  it  belongs  to  the  claaa  of  asymptotic  dlatributiona  of 
tha  extreme  (analleat)  value*  and  nay  adequately  reflect  the  dietribu- 
tion  of  the  proposed  design  pareaMtar,  nenely ,  tha  tine  to  firat  failure 
in  conplicated,  nulticonponant ,  fatigue-critical  • true  turn*.  Moreover, 
tha  Weibull  dietrihutlon  does  reflect  an  increasing  hatard  rata  that 
done  eeen  to  be  a  logical  argument  for  the  fatigue  proceaa.  the  two 
paranetare  conpoaing  tha  atatiatlcal  nodal  represent  fatigue  variability, 
ehape  paranetar  a,  end  the  fat'gua  performance  of  the  etructural 
detail — tha  acale  paranetar  0.  It  waa  decided  that  with  the  aaeunp- 
tion  that  both  paranetar*  are  unknown,  conputed  analyaee  would  be  overly 
conservative  and  a  moro  acceptable  solution  could  raault  if  sons  valua 
for  one  of  the  paranetare  was  established.  Tha  acale  paranetar  was 
known  to  ba  influenced  by  all  tha  variables  of  eatarial,  geometry,  loada, 
environment,  ate. — and  therefore  was  obviously  not  tha  appropriate 
paranetar  for  this  exploratory  task.  Consequently,  effort  was  concen¬ 
trated  on  establishing  the  behavior  of  tha  shape  paranetar  under  actual 
fatigue  conditions.  Unfortunately,  available  fatigue  data  have  been 
almost  entirely  obtained  f roe  testa  with  groups  of  limited  else,  usually 
lovolvlx^  fren  two  to  five  specimens .  Tha  results  of  these  teste  taken 
on  their  own  merits  ere  worthless  for  assessing  tha  typical  parameter 
of  the  parent  popu'mtion,  but  it  wee  assumed  Chat  tha  evaluation  of  a 
large  mass  of  data  obtained  from  numerous  fatigue  tests  would  demonstrate 
the  possibility  of  establishing  tha  shape  paraateter.  This  task  was  of 
sufficient  magnitude  th»t  it  was  decided  to  limit  tha  study  to  only  tha 
aluminum  alloys  and  to  teat  tha  feasibility  of  a  reliability  approach 
on  aluminum  aircraft  structures. 

A  comprehensive  literature  search  of  the  available  aluminum  fatigue 
performance  data  wee  conducted,  a ad  acre  than  2,000  groups  of  data 
represent  lug  approximately  11,000  specimens  were  summarised  and  analysed. 
It  should  ba  noted  that  not  all  tha  groups  were  Independent  of  each 
other  since,  on  a  minority  of  occasions,  a  few  smaller  group#  ware  pooled 
to  form  larger  group#  to  facilitate  evaluation  of  certain  specif led  vari¬ 
ables  within  the  data  sample.  Mora  data  were  collected  during  tha 
literature  search  than  were  actually  analyxad;  same  obvious  trends 
observed  during  the  reduction  of  the  previously  mentioned  groups  of  date 
rendered  unnecessary  an  analysis  of  all  the  data  collected.  These 
observed  trends  concerned  data  from  rotating  banding  taste,  data  on  hand 
forging*,  and  data  from  testa  on  uaaotched  epee  Inane .  The  flret  two 
typea  of  data  both  demonstrated  considerably  greater  scatter  then  was 
noted  for  typical  aircraft  structures,  whereas  the  third  type  showed 
notably  lower  scatter.  Similar  trends  can  be  ot^erved  by  examining  the 
tabulated  results  praeentsd  in  Kef.  12. 
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b.  Analytical  Procadura 

Tha  analytical  technique  adopted  for  aat  lea  ting  the  WaDrall  shape 
par  Mat  at  a  involved  tha  numerical  aolutloe  for  die  aKteM  likeli¬ 
hood  set  law.  te  0CL1) ,  which  waa  cooa&derad  to  provide  either  the  beet1 
estimate  or  nearly  ao  tor  this  toot  and  Indeed  for  ooet  other  applies- 
t Iona .  Tha  sat  Lea ting  aquations  for  both  conplete  end  censored  aeeplea 
given  in  Rat.  13  ara  not  explicit  and  had  to  he  solved  war ically  by 
tha  Mawton  ltarativa  procadura,  which  required  the  uee  of  a  coepoter . 

A  comprehensive  coeputar  prog ran  waa  genarated  that  ami? onpeeeea 
aona  of  tha  many  variables  contained  in  tha  taat  data.  For  example, 
tha  variablaa  of  dif faring  materials,  gaoatatrlee,  loading  coodltlona, 
laboratory  typaa  versus  full-scala  apecineaw,  ate.,  have  been  identified, 
and  tha  computer  output  has  been  arranged  to  provide  this  Inf ores t ion. 

A  breakdown  of  tha  11  digit  description  coda  used  in  the  prog ran  la  given 
In  Appendix  II.  Tha  initial  values  generated  for  each  date  group  were 
point  9vt  lata  tea  of  tha  shape  parameter  1  obtained  f  ran  the  unbiased 
aaxinue  likelihood  eat laates  of  *,  where  oel/tc,  end  the  uabiaeed 
MLE  of  tha  scale  paranater  or  characteristic  life  t.  The  value  of 
characteristic  life  has  bean  Included  to  coeiplete  the  description  of 
the  data  only  and  has  not  been  used  for  any  ensuing  analyses . 

A  study  of  thaaa  point  estimates  of  tha  shape  parameter  showed  that 
several  data  resulted  in  vary  small  values  of  a,  signifying  high  scat¬ 
ter.  A  closer  scrutiny  revealed  e  consistent  trend  in  those  date, 
namely,  the  existence  of  stronger  specimens  with  correspondingly  higher 
eye lee- to- failure  within  what  may  be  called  e  typical  sampling  of 
fatigue  data.  Tha  raaults  studied  would  auggaet  that  thaaa  high  valoae 
occur  within  tha  upper  10S  hand  of  the  failure  times  of  a  fatigue 
sample.  The  Welbull  model  cannot  accept  tha  existence  of  high-time 
outliers  and  when  operating  on  data  containing  this  phenomenon  attempts 
to  account  for  than  by  first  overconpsmeat lag  and  than  indicating  very 
low  tine  failures.  To  Illustrate  this  point  a  few  simulated  groups  of 
data  ware  analysed  and  tha  results  given  In  Table  I. 

An  examination  of  examples  1  and  2  from  Table  I  shows  that  the 
bulk  of  the  data  waa  kapt  constant  but  that,  in  the  firet  case,  one 
specimen  had  a  Ufa  an  order  of  nagnltuda  above  the  average,  whereas 
the  second  example  bed  an  early  failure  an  order  below  the  average. 

It  waa  noteable  that  example  1  with  the  high-time  outlier  showed 
extremely  high  if  cat  tar  with  correspondingly  low  cyclic  values  at  the  10X 
or  less  failure  probabilities,  but  example  2  with  the  low-time  failure 
waa  treated  es  moderately  high  scatter,  with  computed  probabilities  of 
failure  that  ware  aore  consistent  with  the  data. 


She  word  "beat",  as  used  here,  has  a  mathematical  definition.  An 
estimator  la  beat  If  it  p oaeeaeea  minimum  variance  among  all  asti- 
mators.  It  can  be  shown  that  tha  KLX  is  tha  beat  asymptotically 
normal  estimate. 
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TMa  problem  doM  not  disqualify  the  WeiS'll  distribution  aa 
untypical  of  f«tl|M  data;  it  neraly  suggests  that  tha  nloorlty  of  data, 
cantata  lag  apocinana  that  vora  auch  stronger  than  tha  ewer  ago,  should 
ho  co— ored  to  —elude  than  a  apacl— ,  which  obviously  cons  f  ran  s  dif¬ 
fer— t  population,  that  la,  have  a  differ— t  seals  para— tar  I.  This 
process  is  by  so  noons  uapreced— ted ,  as  Walhull  hi— alf  baa  observed 
(14  and  15) :  "So—  of  tha  spocl—  —7  aadurs  a  a— her  of  otr— a  cycles 
auch  larger  than  —  el— ted  f roe  eha  distribution  function"  and  "there 
la  00  a— ae  trying  to  fit  a  singla  function,  nor— 1  or  not,  to  tha  con- 
plats  distribution.  Por  this  ra— on,  tha  data  have  be—  c —sored , 
choosing  a  point  of  truncation  sail  balow  tha  'bn— '  (p  •  90X)." 

It  should  ba  reanph— lead  that  tha  bast  —  tins ea  can—  f roe  a 
consider at  ion  of  all  tha  r— ulta  in  a  data  sat  — copt  for  thoo a  groups 
containing  spocl—  fr—  nixed  populations  uhan  it  would  ba  erromao— 
to  fit  ell  the  date  to  on#  distribution. 

As  a  result  of  tha  previous  argueentc  a  criterion  was  established 
In  tha  HLS  computer  progr—  such  that  whs— ear  high- tine  outliers  wars 
involved,  a  raltarativa  estimation  procedure  was  conducted  to  ranova 
their  influence  on  the  shape  para—  ter.  The  technique  that  was 
developed  is  beat  d— crlbsd  by  raferriig  to  tha  —  plos  given  In 
Table  IX.  The  three  groups  of  a  I— lated  data  pr— tad  are  identical , 
except  that  tha  second  group  cental—  two  additional  high-tine  out liar • 
that  ere  —  order  of  magnitude  above  tha  —jorlty  of  tho  dote,  end  tbs 
third  group  with  two  low-ti—  failures  an  order  of  magnitude  balow  tha 
average  of  tho  renal alng  results.  In  —pis  1,  tha  — cl— ta  of  tha 
reciprocal  shape  perssatsr  fall  within  She  acceptable  range  of  t  <  0, 50 
and  an  —  c— wring  was  at  tamp  tad .  however,  both  examples  2  end  3  con¬ 
tained  data  that  ram ltsd  in  —  tl—  tea  of  the  reciprocal  shape  param¬ 
eter  that  warn  above  tha  arbitrary  threshold  and  that  con— quant ly  ini- 
tinted  tha  ca— orlag  procedure.  It  w—  noted  that  in  axanpla  2  tha 
rejection  of  tho  first  high-el—  outlier  did  not  cause  much  fluctuation 
in  As  — tl— ta  of  *,  but  with  rejection  of  tha  second  outlier  a 
dr/a— tic  chaaga  in  s  w—  affected.  (8—  second  and  third  —  ti—  ta, 
where  tho  ratio  of  second  —  ti—  te/thlrd  — ti— ta  >  1.50).  Further 
censoring  of  tho  data  verified  this  — tl— ta  of  tha  reciprocal  shape 
para— ter  —  stable  aed  therefore  acceptable.  In  ague aple  3,  censoring 
wee  again  attempted  been— a  of  thw  high  scatter  indicated  by  tho  initial 
— ti— ta .  however,  this  scatter  was  tha  rasult  of  tha  two  low-ti— 
failures  and  tha  censoring  of  any  longer  lived  epee  ins—  could  not 
improve  on  eha  initial  — eh— tn,  which  consequently  ra— land  unchanged . 

In  —a ry,  Ala  heuristic  procedure  was  —ad  to  calculate  — wl— 
likelihood-point  —  ti— t—  of  the  data  shape  para— tars,  where — for  the 
large  —jorlty  of  As  groups — ell  As  results  — ra  utilised  in  the 
analysis ,  hue  for  the  ra— lning  minority  of  groups  so—  censoring  of  tha 
data  wao  necessary  to  obtain  reliable  — ti— tea. 

Aa  additional,  — ra  theoretically  acceptabls  — tl— tloc  procedure, 
which  also  circa— tod  Aa  problaa  of  Aa  higb-ti—  outlier,  has  bean 
investigated  and  a  progr—  developed  around  Ala  — tl— tor.  Tha 
a— wing  A— lysis  —  ti—  tea  Aa  Velbull  shape  para— ter  after  considering 


only  the  first-  And  second -ordar  ad  failures  fron  mjt  group  of  dots  of 
any  sit#  n.  Th*  argument  for  thin  eat ins tor  woo,  as  b of or a,  the 
groat  lntoraot  of  the  alrplaiM  doalgnar  ragirflag  th*  an  pec  tad  bUI— 
fatigue  performance  of  a  structure .  Furthermore,  It  wee  known  that  tho 
ext  meal  valuaa  of  a  population  t  and  ad  toward  th*  Ifalbwll  distribution 
in  behavior,  and  tha  axanlnat ion  of  tha  Initial  failuraa  ualag  thla 
a tat  1st  leal  nodal  appaarad  just  If lad.  Thla  r sport  has  uaad  tha  ualquo. 
llnaar,  unhlaaad  invariant  estimator  of  tha  Walhull  shops  parameter  for 
tha  two-ordarad-fallura  aatlantion  procedure,  naaaly t  a  derivative  of 
tha  aaxinun-likalihood  ast  lira  tor  of  r  (5 1/e),  fran  consideration  of 
tha  first  and  sacond  of  n  total  failuraa  [Mann  (16)). 

Ooa  furthat  problem,  naasely ,  how  to  avaluato  tha  estimates  of  tha 
shop*  parameters  raoultlng  fron  tha  analyses  of  tha  aaas  of  eollactad 
fatigua  data,  waa  investigated .  it  oust  ba  r an  nab ar ad  that  tha  vast 
majority  of  tha  data  cans  fron  a  multiplicity  of  tast  groups  of  snail 
sample  also  that  wars  of  llaitad  use  individually  but  that  collactlvaly 
wars  of  aufflclant  auMis  to  auhstantlata  tha  premise  that  tbair  avnraga 
valua  would  convarga  on  tha  truo  value.  Tha  problan  of  weighting 
aatinataa  fron  groups  of  differing  aanpla  aits  waa  wlaad  in  two 
separate  ways.  It  had  baan  noted  that  groups  of  a  aanpla  aits  of  two, 
throe,  four,  and  five  epaciaane  appaarad  to  ba  sufficiently  numerous 
that  a  soparate  traatnant  of  each  aanpla  alas  waa  feasible  and  conse¬ 
quently  it  waa  attempted.  It  waa  later  observed  that  in  Naan's  recant 
paper  (17)  a  procadura  for  weighting  unbiased  aatinataa  fron  groups  of 
dlffarlm  ■ 1b a  had  bean  uaad;  this  procadura  was  therefore  Incorporated 
into  tha  investigation,  thereby  providing  tha  naans  for  generating  an 
unbiased  as tins to  of  tha  walghtad  average  valua  of  tha  shape  parameter 
for  tha  mass  of  fatigua  performance  data.  Ur-fortunately,  it  waa  not 
possible  to  weight  in  tha  ainorlty  of  data  of  email  aanpla  alas  that  had 
been  censored  either  during  tasting  or  by  tha  high- tine- out liar  censor¬ 
ing  process  described  earlier.  This  consent  applies  only  to  tha  MLS 
procadura  and  raaulta  because  two  requisites  of  tha  weighting  process, 
namely,  a  knowledge  of  the  bias  and  variance  values  for  censored  data, 
wore  unknown.  Monte-Carlo  method*  had  base  uaad  to  generate  sons  of  tha 
desired  valuaa,  but  further  raaulta  are  necessary  to  establish  a  compre¬ 
hensive  listing  that  can  than  bo  Incorporated  Into  tha  weighting 
procedure.  This  limitation  in  of  no  real  consequence,  however,  because 
only  a  few  groups  of  data  were  affected  and  unaccounted  for.  Inceptions 
to  this  limitation  occurred  whenever  the  group#  of  censored  data  were 
sufficiently  large  that  asymptotic  theory  could  obviously  ba  applied. 

More  detailed  descriptions  of  tha  analytical  techniques  mentioned 
in  the  preceding  paragraphs  are  given  In  Sac.  IV. 
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2. 


MOULTS  OS  TMORKTIwAL  STUDIES  OH  POPULATION 
An  nalutloi)  of  tho  aanlmum- likelihood  aottaatiag  equations  of 
co— late  or  cortoln  censored  dtu  revealed  chat  tho  following  quantities 
aoro  parameter- free,  itTuitTi,  Independent  of  the  par  ant  population  paras- 
•ttri  (raffr  to  Appendix  I): 

>»»! 


fMKWiOr' 

i  (refer  ti 
U  •  |/o  t 

v  -  <{/•>; 
w  -  (♦/•)• 


Tha  unknown  population  parameters  a  and  8  ara  the  true  Walbull  a  ha  pa 
and  acala  values  respectively ,  wharaea  8,  ft,  and  i  ara  tha  corre¬ 
sponding  MLS ,  with  ft  tha  valua  whan  a  in  bom  and  I  whw  a  la 
unknown. 


because  tha  distribution  of  tha  atatlatlc  If  ia  known  exactly  from 
theory ,  but  thoao  of  U  and  V  vara  unknown,  computer  sinulatlon  tech¬ 
niques  wara  uaad  in  conjunction  with  MLS  procedures  to  obtain  anplrlcal 
d 1st rlbut Iona  fov  tha  unknown  quantities  aantlonad  above ,  Sana  maser  1- 
cal  results  ara  (Ivan  In  Tab las  III,  IV,  and  V;  thaaa  results  were 
plottsd  on  normal  probability  papar  (figs.  4,  5,  and  6)  to  lndlcata, 
graphically,  tha  davlation  of  tha  important  statistics  1/U,  V,  sad  W, 
from  the  "par fact"  valua  unity  for  various  groups  of  data  of  finite 
aanpla  else. 

Figure  A  shows  the  reciprocal  transformation  of  tha  aaplrlc  distri¬ 
bution  of  ft/a.  This  reciprocal  of  U  was  consider ad  b acauaa  it  had 
racalvad  recant  attantlon  la  tha  litaratura  (18).  Tha  distribution  1/U 
for  tho  sample  aiaa  n  •  2  was  also  known  exactly  from  thaory  end  has 
baaa  included  tc  demonstrate  tha  "goodness"  of  tha  sinulatlon  rooulto 
obtained  in  this  study. 

Figure  S  shows  tha  empirical  distribution  of  tho  NUB  of  tha  Welbull 
seals  parameter  0,  without  tho  booeflt  of  any  prior  knowledge  regard¬ 
ing  tha  shape  parameter  a.  It  is  obvious  that  >nder  this  raatrictlon 
a  designer,  given  tha  raaults  of  a  teat  on  flvo  or  fewer  spec inane ,  nuet 
neks  sn  unfortunate  choice  between  low  confidence  In  average  or  charac¬ 
teristic  life  prediction  or  emcees  weight.  Figure  6  shows,  however, 
that  given  a  prior  knowledge  of  tho  smact  valua  of  o,  samp  lea  of 
only  aae  nay  bo  adequate  for  obtaining  beth  high  confidence  levels  sod 
reasonable  weight  simultaneously.  Some  stags  between  thaaa  aatramaa  of 
no  knowledge  and  complete  knowledge  of  tho  aha pa  parameter  would  be 
representative  of  reality. 

Tab las  Ill  and  IV  include  roaults  to  provide  eons  insight  into 
high-time,  single-stage,  type- II  censoring .  Type  II  censoring  involves 
tha  cessation  of  tooting  after  the  occurrence  of  a  predetermined  nuaber 
of  failures,  whereas  type  I  censoring  refers  to  that  fora  of  tastii*  In 
which  the  rnnlra  number  of  load  applications  is  pradeteimlnad  and  thus 
results  in  the  cessation  of  testing  either  by  epacimea  failure  or  by 
the  aspiration  of  tho  preoat  tins  period.  Single-stage,  type  II  censor¬ 
ing  wan  of  paramount  Interest  not  only  because  it  typified  the  nsjorlty 
of  the  censored  data  groups  but  also  because  It  was  tha  nodal  for  tho 
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two-orderad-fallura  estimation  procedure  used  in  this  r apart .  Appendix  I 
aatabllahaa  that  under  thla  fora  of  censoring  1  A' ,  V,  and  V  ara  parameter  - 
fraa. 


Tab! an  III  and  IV  give  'hr  distributions  of  i/U  and  V 

for  two  exaaplae  ,  "mealy,  tha  complete  caaa  of  thraa  falluraa  trim  a 
amp  la  of  thraa,  and  tha  canaorad  caaa  rf  thraa  falluraa  f  ron  a  cample 
of  five  ■ pat In ana .  It  can  ba  aaan  that  tha  acala  paranatar  •  la  aora 
accurately  aatlmatad  from  tha  canaorad  caaa  of  tha  thraa  weakest  In  five, 
but  tha  ahapa  pa  rant#  tar  a  la  aatlmatad  aqually  wall  from  tha  ana  liar , 
comp lata  aaapla.  Now  conaidar  tha  example  whar  all  five  failura  tlmaa , 

In  tha  aaapla  of  five,  war#  known,  rathar  than  only  tha  thraa  waakaat  In 
flva.  Tha  cables  show  that  tha  aatlnatlon  of  a  la  laprowed  by  tha 
incraaaad  knowledge,  but  that  tha  aatlmata  of  0  la  lnprovnd  only 
slightly.  Thasa  11m It ad  exaaplae  do  raflact  soma  Important  traada,  auch 
aa  tha  improvement  in  aatlnatlon  of  acattar,  or  a,  gainmd  by  tba  know¬ 
ladga  of  both  low  and  high  tlmaa  to  failura  rathar  than  juat  ona  of 
tha  extrema*.  Furthermore,  It  ammaa  that  tha  third  failura  In  a  aampla 
of  fiva  la  a  battar  indicator  of  central  tendency,  or  of  0,  than,  for 
axanpla,  la  tha  eecond  failure  in  a  aaapla  of  thraa. 

Figure  7  la  a  Joint  presentation  of  computational  raaulta  obtained 
during  thla  study  and  of  raaulta  given  in  a  recant  paper  (10).  Snap ling 
error  In  tha  maxima*- likelihood  aatlmata  of  a  was  neaaurad  by  tha  plot 
of  tha  variance  and  of  tha  expected  loaa  of  1/U  (l.a.  «/ft)  varaua 

aampla  alia  (complete  samp lea ) .  Soma  additional  (rectangular)  polnta 
were  located  on  this  plot  end  give  she  expected  loss  of  «/*.  The  letter 
estimate,  called  the  simplified  method  of  aomeota,  la  obviously  more 
easily  obtained  than  ft,  but  la  also,  as  a  hewn  In  Fig.  7,  leas 
accurate.  In  fact,  Raf.  10,  which  review*  all  estimating  procedures  for 
tha  Welbull  parameters,  obtained  raaulta  showing  that  in  moat  caaaa  tha 
Maximum- likelihood  estimates  had  the  least  expected  loss  and  hence  the 
noet  accuracy  of  all  reviewed  estimators. 

3 .  RESULTS  OF  DATA  ANALYSES 

a.  General  Observations 

Appendix  II  itanlaaa  tha  individual  characteristics  of  each  of  tha 
collected  data  groups.  Data  for  2,000  groups  representing  11,000  epeci- 
mm  are  Hated  and  include  fatigue  performance  data  on  tha  well-known 
aluaimm  alloys  of  2024,  7075,  and  7178,  together  with  soma  of  tha 
leaser  used  alloy*  such  aa  2014  and  a  variety  of  tha  European  alloys 
Ilka  DTO687A,  DTD610B,  ate.  Specimens  formed  from  sheet,  plate,  bar, 
ax true Ion,  and  forging  ara  rapr as anted  in  several  configurations  rang¬ 
ing  from  tha  a  lapis  uncotchad  and  notched  types  to  structural  simulator* — 
such  a a  mechanically  fastened,  or  bonded,  or  *pot-welded  joints — to 
structural  component*  and  full-scale  structures.  Tasting  conditions 
encompassed  axial  loading,  both  tana lie  and  compressive;  flexural  and 
rotating  banding;  both  constant-amplitude  and  variable-amplitude  tasting 
techniques ;  and  ranged  from  lew  cycle  (below  10?  cycles)  through  to  high 
cycle  (above  10*  cycles)  fatigue  performance  tasting.  It  is  also  notable 
that  tha  great  majority  of  tha  data  comae  from  groups  of  a  sample 
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■lie  of  two  to  five  (fwctatn*  and  that  only  a  ana 11  alacrity  of  groups 
contain  20  or  aore  specimens. 

A  stud"  of  tha  raaultant  point  aetlmtes  of  tha  a  ha  pa  parameter 
given  in  Appendix  II  revealed  soma  obvious  trooda.  Thaaa  concerned  tha 
large  acattar  demons tratad  by  rotating  banding  taata  and  by  fatlgua  data 
on  hand  forgings,  with  both  groups  raaponalbla  for  earns  of  tha  largaat 
acattar  valuaa  recorded  (a.g.  I  <  0.3).  It  waa  also  no tad  that  several 
caata  on  axially  load  ad ,  adge-ootchad  specimens  demonstrated  consistently 
lowar  acattar  than  nornal.  Further  investigation  revealed  that  tha 
speclmmne  uaad  In  thaaa  taata  had  baan  manufactured  with  exceptional 
care  f row  a  controllad  stock  of  notarial  to  as  to  ainlwlsa  tha  offsets 
of  fatlgua  acattar  on  tha  anculng  taata  (Haifa.  19,  20,  and  21). 

Tha  obvious  c rands  just  nantlonad  dawonatratad  tha  existence  of 
aana  data  within  tha  col lac tad  aass  of  fatlgua  raaulta  that  could  ba 
dlaquallflad  as  untypical  of  tha  dlaparalon  expected  f row  aluninua 
airplane  atructuroa.  It  also  suggest ad  tha  advisability  of  a  cloaar 
scrutiny  of  all  tha  data,  which  waa  subsequently  undartakan. 

b.  Factors  Influencing  Fatlgua  Scattar 

Tha  tvo-ordsr-statlstlc  aa tins tor  in  conjunction  with  a  weighting 
technique,  which  ara  both  axplainad  in  Sac.  IV,  waa  uaad  to  axawlna 
savaral  of  tha  known  varlabloa  wiehle  tha  fatlgua  parfomanca  data.  Tha 
raaulta  of  tha  analysas  hava  baan  tabu la tad  in  Tabla  VI  and  plottad  In 
Figs .  8  thrsajh  14,  which  reveal  several  ietaraating  aspects  of  tha 
fatigue  parfomanca  of  aluQlmaa.  It  is  apparent  free  Tabla  VI  that  tha 
ltowixad  varlabloa  ara  not  wholly  independent,  aa  any  one  tast  result 
way  ba  Included  In  several  of  tha  11a tad  catagorlaa.  Nevertheless,  ha 
comparisons  nada  ara  justified  and  demons trata  that  with  few  exceptions 
fatigue  data  cows  f raw  a  parent  population  having  a  Walbull  a ha pa  param- 
atar  In  tha  rang*  *  <  a  <  5. 

Tha  Initial  run  comisted  of  all  tha  gathsrsd  data,  which  contalnad 
data  that  had  vary  high  acattar  (a.g.  hand  forglnga,  uanotchad  rotating 
banding,  otc.)  and  othar  data  with  vary  low  acattar  fro*  carafully 
controllad  taata.  Tha  oppoolta  trends  of  thaaa  data  t and ad  to  ba  salf- 
cancellqg ,  and  tha  raaultant  waightad  svaraga  aatlnata  of  k  (l.a.  *") 

waa  a  fairly  sc cap tab la  valua.  Tha  nest  rasult  * - — 4  — 1  only  that  data 

suapactad  of  being  at  sows  diaagraawant  with  tha  ganaral  trends,  that  la, 
tha  it ana  Just  enumerated.  Thaaa  raaulta  wore  plottad  in  Fig.  8  and 
describe  tha  grantor  scatter  that  waa  axpactad.  Tha  third  result , 
represented  in  Tabla  VI  and  plottad  in  Fig.  9,  shows  th*  affect  of 
isolating  only  that  data  found  to  typify  tha  fatlgua  acattar  experienced 
by  aircraft  structural  components.  Thaaa  data  included  raaulta  frow: 

1)  All  tha  Invoatigatad  aluwinua  alloys 

2)  All  tha  notched  configurations ,  ranging  frui  tha  siwple  wooo- 

llthlc  notched  specimen  to  tha  full-scale  structural  component 

3)  both  axially  loaded  and  f laurel ly  loaded  taste 

4)  loth  constant-  and  variable-amplitude  taata 

5)  Tha  fatlgua  performance  range  of  10?  to  10*  cycles 
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It  Is  obvious  chat  the  criteria  lieted  above  cover  a  broad  rang*  of 
data,  and  the  subsequent  result e  and  figures  are  presented  to  eubetan- 
tlate  this  conclusion.  Before  proceeding,  however,  It  should  be  noted 
that  with  the  rejection  of  specific  data  that  had  demonstrated  untypicaUy 
large  acattar  (e-g.  notched  rotating  bending  testa,  results  from  hand 
forgings,  low-amplitude/hlgh-llf a  fatigue  dsta,  compression-compression 
test  results,  end  dsta  on  bonded  joints,  together  with  data  ehowlt* 
untypically  low  scatter  (e.g.  fatigue  teats  on  unnotched  epee  lamia  and 
also  on e  particular  set  of  notched  data  that  had  bean  extremely  care¬ 
fully  machined  to  specifically  alnimlte  scatter]),  the  weighted  average 
estimate  of  <  la  reduced  from  *'  •  0.26  for  all  data  to  *'  <•  0.22 
for  the  apecifled  data. 

Figure  10  coepai  r*  the  ciaaulative  frequencies  of  th*  three  alumi¬ 
num  alloys:  2024,  707b,  and  7178.  Only  thoae  data  meeting  Che 
apecifled  criteria  (listed  in  the  preceding  paragraph)  have  been  analysed, 
end  it  can  be  seen  that  there  are  no  obvious  differences  In  the  distri¬ 
butions  of  these  material*.  The  7178  curve  ahowa  some  deviation  over 
pert  of  the  range,  but  It  should  be  noted  that  the  sample  alas  was  com¬ 
paratively  .limited.  Figure  11  compares  the  differing  specimen  configura¬ 
tions  of  simple  monolithic  notched  email  structural  simulators,  such  as 
lap  joints,  and  the  full-scale  structures  and  structural  components . 

Again  no  real  differences  were  noted  between  the  various  geometries, 
and  a  comparison  with  the  data  in  Fig.  10  showed  that  the  cumulative 
frequency  plots  were  similar  for  the  differing  variables  of  material 
and  geometry.  Figure  12  compares  the  results  of  constant-amplitude 
and  variable-amplitude  testing  techniques,  and  it  can  be  seen  that 
verieble-amplltuda  tasting  does  ahow  slightly  lower  scatter.  It  should 
be  noted  that  the  sample  of  variable-amplitude  data  contained  considerably 
fewer  groups  than  the  constant-amplitude  data,  but  this  alone  might  not 
account  for  all  the  difference.  A  closer  study  of  the  variable- 
amplitude  data  revealed  that  one  set  of  results  obtained  from  Ref.  22, 
showed  behavior  different  from  the  observed  trend.  It  had  been  noted 
during  the  course  of  this  study  that  the  ratio  of  time  to  aacond  failura 
versus  tins  to  firs':  failure  was,  on  average,  approximately  1.4 
[i.a.  E(X2/Xi)  :  1.4],  but  for  these  data  it  was  noted  that  tha  initial 
thrae  or  four  failures  were  all  within  that  ratio  and  that  tha  first 
and  aacond  failuraa  were  right  on  top  of  each  other.  An  estimate  of 
this  sample  of  data  gave  a  value  of  R’  *  0.12,  which  was  considerably 
below  the  rest  of  the  variable-amplitude  data,  where  »  0.19.  How¬ 
ever,  because  no  obviously  acceptable  reason  was  available  for  the  rejec¬ 
tion  of  these  data,  they  were  retained  In  the  variable-amplitude  sample. 

Figure  13  plots  the  simple  monolithic  unnotched  data,  bonded  joints 
date,  end  low- amplitude /high-life  fatigue  performance  date.  It  is 
obvious  that  these  raaults  are  substantially  different  from  those  given 
by  tha  qualified  data,  and  they  ware  consequently  disqualified  as  being 
untypical  of  the  variability  expected  from  alunlntm  aircraft  structures. 
Furthermore,  it  could  be  argued  that  laboratory- type  unnotched  specimens 
were  sufficiently  different  from  structural  details  containing  a  multi¬ 
tude  of  holes  and  discontinuities  that  their  worth  as  a  test  specimen 
wee  doubtful.  Tha  bonded  data  came  fro*  same  older  report! ,  which 
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exhibit ad  iom  inconsistencies  in  the  bonding  process  end  thus  contained 
test  results  with  extremely  early  failures  and  therefore  large  scatter. 
Finally,  low-amp 11 tude/hlgh-cycle  data  reflect  the  effects  of  the 
endurance  Unit  where  the  shallow  elope  of  die  S-H  curve  tende  to 
exaggerate  any  differences  in  the  test  loading.  The  airplane  atructure 
is  designed  for  sore  finite  lifetimes,  and  the  exaggerated  scatter  of 
the  very-long-lived  teats  would  be  unrepresentative  of  this  design 
parameter .  The  cowmen ts  enumerated  above  add  sons  further  justifica¬ 
tion  for  neglecting  these  data  froa  the  estimation  of  the  typical  shape 
parameter . 

Having  observed  the  large  scatter  of  the  data  from  high-life  testing, 
the  remaining  mass  of  data  wcj;  examined  to  check  for  any  trends  show¬ 
ing  increase  in  scatter  with  increase  in  cyclic  lives.  The  range  of 
data  between  70  and  10 6  cycles  of  constant-amplitude  loading  was  broken 
into  five  ranges  that  were  considered  of  general  Interest  to  the 
designer,  and  estimates  of  the  shape  parameter  were  conducted  for  these 
ranges.  Figure  1A  shows  these  results  and  demonstrates  that  within  the 
broad  range  of  70  to  10°  cycles  fatigue  scatter  does  not  increase  with 
cyclic  life.  However,  above  and  below  this  range  scatter  was  noted  to 
increase,  which  could  be  expected  considering  the  fist  tails  at  alther 
end  of  the  S-N  curve. 

These  analyses  of  the  mass  of  collected  data  have  revealed  several 
important  facts  regarding  the  fatigue  scatter  of  aluminum  alloys.  They 
have  shown  that,  with  certain  specific  restrictions,  fatigue  scatter 
did  not  vary  with: 

1)  The  material  of  the  aluminum  alloy  specimen 

2)  The  geometry  of  the  notched  specimen 

3)  The  cyclic  life  of  the  specimen 
A)  The  fatigue  test  machine 

5)  The  type  of  loading  (whether  constant  or  varlabla  amplitude) 

Actually,  a  trend  toward  marginally  lower  scatter  was  observed  froa 
variable-amplitude  data  when  compared  with  constant-amplitude  data, 
but,  because  the  variable-amplitude  data  contained  considerably  fewer 
samples,  further  work  would  be  necessary  to  make  any  more  of  this 
slight  trend. 

It  is  realized  that  these  results  do  not  substantiate  some  published 
conclusions  (Refs.  23  and  2A)  that  were  baaed  on  a  limited  amount  of 
test  date  of  which  many  were  judged  to  be  untypical  of  aluminum  struc¬ 
tures.  It  should  be  emphasized  that  conclusions  mads  about  population 
parameters  oust  be  generated  from  a  background  that  is  sufficiently 
large  to  reduce  the  problem  of  both  systematic  and  sampling  error,  and 
it  Is  believed  that  this  current  investigation  has  taken  a  major  step 
toward  this  goal. 

c.  Comparison  of  Empirical  and  Theoretical  Sampling  Distributions 

During  the  analysis  of  the  data  an  observation  had  been  made  of 
the  great  preponderance  of  samples  of  size  two,  three,  four,  and  five 
specimens.  These  Individual  groups  were  examined  more  closely  to  cir¬ 
cumvent  the  problem  of  weighting  information  from  groups  of  differing 
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sample  size.  Table  VII  presents  the  results  of  the  three  analyses 
conducted:  (1)  the  MLE  procedure  for  data  that  Include  high- time- 

outliers  ,  (2)  the  ad  hoc  high-tlme-out lier  censoring  procedure,  and 

(3)  the  two-ordered-failure  estimator.  The  tabulated  results  demonstrate 
the  improvements  in  the  estimates  of  *■  that  were  obtained  by  the 
censoring  of  data  containing  high-time  outliers.  The  samples  of  size 
two  could  not  be  censored,  fcr  obvious  reasons,  and  therefore  do  not 
alter,  but  the  remaining  groups  all  reflact  lwer  ecattar  with  rejection 
of  the  high-time  outlier  and  follow  a  consistent  trend  regardless  of 
the  censoring  method  used . 

Figures  15  through  18  show  plots  of  the  percentage  frequency  end 
cumulative  frequency  of  2  for  samples  of  two  to  live  specimens.  Those 
empirical  distributions  have  been  compared  with  the  theoretical  cumu¬ 
lative  sampling  distributions  that  were  determined  from  Eq.  (TV-13)  by 
using  the  value  of  £'  for  the  1,298  groups  of  qualified  deta,  that  is, 

-  0.224.  (See  Table  VII.)  All  four  figures  exhibit  the  same  trends, 
a  similarity  of  the  empirical  and  theoretics!  curve* ,  with  the  date 
showing  slightly  less  conservatism  than  that  predicted  by  theory,  except 
at  the  higher  percentiles  where  a  croesover  of  ti»e  curves  is  observable. 

To  compensate  for  this  latter  behavior  it  was  necessary  to  adopt 
some  interval  estimate  of  the  shape  parameter,  and  a  951  confidence 
level  was  assumed  to  provide  as  severe  an  interval  as  desirable.  The 
qualified  data  were  of  such  a  considerable  mass  that  the  resultant 
deviation  of  the  collective  estimator  was  very  small,  being  in  tha  order 
of  0.006.  As  nothing  further  was  known  about  the  distribution  of  the 
estimate,  the  Tchebycheff  inequality,  which  leads  to  the  statement 
p  (k  <  t'y)  >  y ,  was  used  for  establishing  the  desired  upper  bound,  i'y. 
This  resulted  in  a  value  of  £'o.95  "  0.25,  and  curves  generated  from 
this  number  have  also  been  plotted  in  Figs.  15  through  18.  These  plotted 
curves  help  demonstrate  the  expected  improvement  in  reliability  to  be 
gainec  by  using  the  interval  value  of  £! '  ,  and  it  can  be  seer,  from 
the  proximity  of  the  theory  curves  of  the  weighted  mean  value  and  the 
internal  value  of  that  the  added  conservatism  given  by  the  interval 

estimate  should  not  result  in  an  excessive  increase  in  the  overall 
"scatfer  factor." 

Figures  19  through  24  are  plots  of  the  percentage  frequency  and 
cumulative  frequency  of  £  for  specific  sample  sizes  from  the  total 
mass  of  collected  date.  These  have  been  included  to  facilitate  compari¬ 
sons  between  the  total  collected  data  and  that  finally  qualified  aa 
typical  of  aluminum  aircraft  structures.  A  comparison  of  these  graphs 
with  Figs.  15  through  18  demonstrate  the  greater  divergence  between  the 
empirical  and  theoretical  curves  for  the  total  data.  This  is  further 
accentuated  by  observing  the  much  closer  correlation  of  the  sample  mean 
values  with  the  average  mean  value  given  in  the  plots  of  the  qualified 
data  in  Figs.  15  through  18.  These  trends  are  further  verification  of 
the  improvement  (i.e.  reduction)  in  systematic  error  that  was  obtained 
through  rejection  of  specific  data. 
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Table  VIII  has  been  included  to  substantiate  the  improvement  gained 
by  the  editing  of  the  tor.al  naan  of  collected  date  to  ait  elude  thoee 
considered  to  be  untypical  or  unrapreaentat lve  of  alumlmst  structural 
behavior.  The  reeulte  tabulated  include  the  aaaple  else,  the  number  of 
group*,  the  estimated  value  the  variance  of  the  observed  distri¬ 
bution  of  *;  and  then  two  theoretical  variances,  tha  first  beast!  on 
the  individual  estimate  of  «  and  tha  second  calculated  from  tha  overall 
welgnted  average  eetlmate  *’  obtained  from  all  the  qualified  fatigue 
performance  data.  The  initial  results  presented  are  comparisons  of  tha 
observed  and  theoretical  values  for  samples  of  site  two,  three,  four, 
snd  five  obtained  from  an  examination  of  all  the  collected  data.  It  It. 
obvious  from  these  reeults  that  th«>  observed  variance  was  greater  th.o 
that  predicted  by  theory.  The  rMMlnng  result#  were  obtained  from  only 
the  data  that  vara  qualified  aa  typifying  currant  structural  performance, 
and  it  can  be  seen  that  tha  observed  variance  has  now  reduced  to  almost 
the  equivalent  value  as  was  obtained  theoretically .  A  further  compari¬ 
son  of  tha  empirical  variance  with  the  theoretical  variance  obtained 
from  the  weighted  average  value  of  also  shows  s  good  correlation  for 
the  qualified  data  but  a  rather  poor  agreement  for  the  total  collected 
data. 

4.  DETERMINATION  OF  SCATTER  FACTORS 

Table  VII  gives  the  results  of  tha  analyses  conducted  on  fatigue 
performance  data.  The  two  methods  used  for  compensating  for  the 
influence  of  high-time  outliers  within  certain  groups  of  data  resulted 
in  two  slightly  different  astlautes  of  r,  namely:  r  -  0.2416  for 
the  MLE  method  using  all  tha  data  points  and  <’  *  0.2244  for  the  two- 
ordered  failure  estimating  method.  Naturally  these  two  estimators  are 
weighted  differently;  the  deviation  of  the  former  was  approximately 
0.003,  whereas  the  latter  vaa  approximately  0.006.  Using  thsse  values. 
Interval  estimates  of  <  were  established  at  the  93Z  confidence  level 
through  the  Tchebycheff  inequality,  and  the  resultant  values  ware  found 
to  be  almost  identical  at:  o  95  "  0.2375  end  < 95  “  0.2525 

Aa  a  consequence  of  thia  equality  in  the  results  it  was  decided 
to  adopt  the  f lrst-tvo-ordered-failuras  estimating  procedure  for  the 
following  reeaons : 

a)  Ease  of  computation  (does  not  require  iterative  techniques) 

b)  Handles  the  problem  of  the  high-time  outlier  without  resorting  to 

heuristic  techniques 

c)  la  an  order-statistic  approach  to  the  problem 

Having  established  an  upper  bound  for  Che  Welbull  shape  parameter 

at  a'Q  95  -  4.0  (i.e.  -  0.25)  for  all  typical  aluminum  airpLane 

structure,  it  is  possible  to  study  the  effects  of  the  variables  that 
combine  *0  produce  the  "scatter  factor."  In  this  report,  scatter 
factor  is  defined  as  the  ratio  of  the  tested  or  calculated  average 
or  characteristic  life  versus  the  design  or  certifiable  life.  The 
scatter-factor  variables  proposed  in  this  report  are  defined  aa: 
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a)  The  sraple  slae  factor  Sn^,  which  regulates  for  tha  improvement 

In  the  estimates  of  tha  acala  or  location  pa  ran# tar  that  raault 
f row  increasing  tha  taat  Input  information  front  a  single  failed 
spec  la an  to  an  infinite  maaber  of  data  point*. 

b)  The  reliability  factor  S|,  which  providaa  tha  ralatlooahlp 

between  some  level  of  performance  (e.g.  ataan  or  characteristic  Ufa) 
and  tha  performance  at  tha  dee  1 red  level  of  safety  foT  a  fleet  of 
one . 

c)  The  fleet  size  factor  Sjjj,  which  compensates  for  tha  decrees 

reliability  caused  by  varying  the  fleet  exposure  from  a  single 
to  a  multiplicity  of  unite. 

The  proposed  scatter  factor  is  the  product  of  these  three  variables: 

Scatter  Factor  (S)  •  S  *  Ss  •  S„ 

nf  *  *Ht 

The  arguaents  for  this  approach  have  already  been  presented  In 
Sec.  II,  end  the  analytical  definition  of  this  reliability  plan  la 
given  in  Sec.  V. 

Table  IX  has  been  presented  to  illustrate  the  relative  influence  of 
any  one  variable  on  the  overall  scatter  factor.  The  char acterle tic  life 
is  the  Welbull  scale  parameter  g  ,  which  is  estimated  by  using  the 
MLE  equations  on  the  results  of  full-scale  structural  fatigue  taste. 

The  values  presented  for  the  variation  of  sample  else  (i.a.  So.)  include 
the  factor  for  converting  tha  point  estimate  of  p  to  a  one-sided 
lower  bound  estimate  at  the  95X  confidence  level.  It  is  readily 
apparent  from  the  table  chat  only  a  mall  penalty  la  incurred  by  testing 
the  minimum  rumba r  of  full-ecale  developmental  teste.  However,  it  is 
equally  obvious  that  the  reliability  specification  and  fleet  size  ere 
critical  components  on  the  final  aagn~tude  of  the  scatter  factor  and 
therefore  require  careful  consideration  when  establishing  a  daalrad 
performance  goal. 

5.  APPLICATION  OF  RELIABILITY  ANALYSIS  TO  AIRPLANE  STRUCTURES 
a.  Example  Problems 

The  proposed  reliability  analysis  plan  has  been  applied  to  a 
few  structural  details  for  which  both  full-ecale  laboratory  taat 
results  and  actual  service  performance  data  ware  available.  The  air¬ 
plane  considered  was  a  military,  jet-eoglnmd  tanker /transport- type 
aircraft  of  approximately  300, 000- lb  maximum  gross  weight. 

The  first  example  concerned  a  traillng-edge  discontinuity  detail 
on  the  wing  lower  surface.  From  full-scale  structural  taat: 
e  Time  to  detected  crack  -  3,697  spectra  on  starboard  wing  only 
•  Equivalent  flight  time  per  load  spectrum  application  ■  5  hours 
e  Time  to  first  detected  crack  -  18,500  hours 
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A  retrofit  package  planned  to  extend  the  service  life  had  been  Initiated 
end,  aa  of  itnusry  1969,  268  airplanes  of  the  soma  particular  type  had 
been  inspected  and  nod  if  led.  Low  tine  for  a  crack  out  of  this  r.roup  of 
inspected  airplanes  "  2,720  hours. 


A  comparison  of  cha  current  fatigue  evaluation  method  and  the  pro¬ 
posed  reliability  analysis  eathod  is  itemised  bale*#: 

1)  Arbitrary  scatter  factor  of  4 

a  Test  life  -  18,300  hours 

e  Safa  life  -  18,500/4  -4,600  hours 

2)  Reliability  analysis 

e  Test  life  on  starboard  wing  -  18,500  hour?  and  the  port  wing 
>  18,500  hours 

e  Teat  characteristic  life  -  22,000  houra 


from  d  - 


+  (n-nf)  Y 


a 

<nf)J 


1/a 


(V-2) 
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and  from  Tab  la  IX: 


Teet  sample  size  -  1  failure 

sn# 

-  1 

.32 

Float  site  (268  airplanes )  -  536  wings 

Sn 

-  4 

.81 

Reliability  of  mean  Ufa  of  waakaat 

Sfi 

-  1 

1C  2 

75X  reliability  of  waakaat 

SR 

-  1, 

.365 

Mean  life  of  weaken t  of  fleet,  therefore,  -  22,000/7.0 

•  3,140  hours 

75X  reliability  of  weakest  of  fleet,  therefore,-  22,000/8.7 

-  2,530  hours 


However,  severe!  other  structural  details  ware  also  Inspected  end 
modified  during  installation  of  the  previously  mentioned  retrofit  pack¬ 
age,  end  these  results — together  with  that  just  given  in  the  example — 
are  tabulated  in  Table  XI.  The  expected  performance  given  by  the  relia¬ 
bility  analyses  and  based  on  test  livss  of  the  detsila  or  calculated 
performance  based  on  s  test  result,  all  at  the  95Z  confidence  level,  are 
also  compared.  The  several  differing  versions  of  the  reference  tanker/ 
transport  airplane  were  examined  end  are  also  Included  in  the  table. 

It  is  Immediately  obvious  that  the  proposed  scatter  factors  vary  widely 
on  both  sides  of  the  current  value  of  4,  depending  on  the  size  of  the 
exposed  fleet  end,  consequently,  estimate  safe  lives  that  are  different 
from  those  currently  computed.  It  la  also  obvloua  that  the  fatigue 
Improvement  packages  were  installed  at  times  that  were  baaed  on  scatter 
factors  >  4. 


A  second  example  again  concerns  wing  structure  in  which  skin 
cracks  were  discovered  at  a  critical  skln/doublcr  detail  From  the 
full-scala  structural  teat: 

a  Time  to  detected  crack  *  1,498  apectra  o;  etarbe  ing  only 

a  Equivalent  flight  time  per  load  spectrum  applies -'.u:  -  5  hours 
•  Time  to  first  detected  crack  -  7,500  hours 

Low- time  failure  in  fleet  -  655  hours 
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A  cmpirlion  of  the  currant  fatigue  evaluation  aethod  and  cha 
proposed  reliability  analyala  aathod  la  itemised  below: 

1)  Arbitrary  acattar  factor  of  A 

a  Teat  Ufa  -  7,500  hours 

e  Safe  life  -  7,500/4  -  1,875  hour* 

2)  Reliability  analyala 

a  Teat  life  on  starboard  wing  -  7  ,500  houra  and  the  port  wing 
>  7 , 500  houra 

a  Teat  characteristic  life  -  8,920  houra  {(fron  Kq .  (V-2))) 
and  fro*  Table  IX: 

Teat  sample  size  ”  1  failure  -  1.32 

Fleet  size  (500  airplanes)"  1,000  wings  Sg  •  5.63 

Reliability  of  aean  life  of  weakest  S|  *  1.102 

75X  reliability  of  weakest  Sf  »  1.365 

a  Mean  life  of  weakest  in  fleet,  therefore,  •  8,920/8.2 

»  1.100  hours 

a  7SX  reliability  of  weakest  of  f leet , therefore, •  8,920/16.1 

-  880  houra 

It  is  readily  apparent  that  the  order-statistic  application  of 
reliability  analysis  forecasts  a  elninue  service  life  that  la  recog¬ 
nizably  less  than  the  conventional  approach.  Moat  noteworthy  la  the 
closer  approximation  of  actual  service  data  by  this  proposed  approach. 
An  obvious  hazard  of  the  proposed  reliability  analysis  plan  la  its 
total  dependence  on  the  test  result,  with  its  aseueed  equality  of  test 
and  service  environment .  It  is  imperative,  therefore,  that  as  accurate 
and  reliable  an  approximation  to  service  conditions  ss  is  practical 
be  considered  in  establishing  the  characteristic  life  of  the  structure. 
At  this  time,  however,  the  actual  correlation  between  service  and  pre¬ 
dict  ion  seems  promising. 

One  further  point  regarding  the  preceding  examples  concerns  the 
choice  of  degree  of  reliability  of  the  weakest  embers.  If  it  is 
extremely  important  that  the  weakest  members  be  guarded  against  the 
likelihood  of  failure,  some  reliability  interval  other  than  median  or 
mean  should  be  considered.  A  75X  level  of  reliability  was  also  assumed 
ir.  the  preceding  examples  for  demonstrative  reasons  only  and  not  aa  a 
suggested  value  for  the  reliability  interval.  It  can  be  seen,  however, 
chat  even  this  fairly  low  specification  would  have  proved  adequate  for 
predicting  the  failure  in  the  first  example  and  approached  the  perform¬ 
ance  level  of  the  weakest  in  the  second  eomaple  It  should  be  added 
that  the  second  example,  although  a  fatigue  failure,  had  bean  aggravated 
by  poor  initial  assembly,  thus  causing  an  extremely  early  failure. 

b.  Possible  Design  Penalty  Ftcm  Application  of  Reliability  Analysis 

Plan 

Some  additional  work  was  completed  to  provide  a  Measure  of  the 
penalty  that  eight  be  Incurred  by  the  introduction  of  this  proposed 
reliability  analysis.  Fatigue  performance  analyses  were  conducted  at 
several  fatigue  critical  locations  on  the  jet-engined  military  tanker/ 
transport  aircraft  referred  to  earlier  in  tbe  discussion.  A  typical 
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6-hour  ■lesion  profile  was  used  in  conjunction  with  gust  and  maneuver 
data  chat  had  bssn  reduced  from  SAC  VHG  date  a  assured  on  the  reference 
airplane  at  Castle  and  Walker  Air  Force  Bases  (kef.  25).  Figures  25, 

26,  and  27  describe  the  mission  profile,  gust  dr.ta,  and  maneuver  data 
respectively .  Reference  locations  were  selected  free  each  of  the  eajor 
components  of  the  airframe  and  comprehensive  fatigue  damage  analyses 
conducted  to  determine  their  respective  performances .  These  analyses 
were  based  on  the  Faiagren-mner  fatigue  damage  hypothesis  together  with 
stress  versus  cyclic  life  date  obtained  from  fatigue  teste  on  full-scale 
structural  components.  An  example  family  of  S-N  curves  le  shown  in 
Fig.  29.  Figure  29  was  plotted  from  the  results  of  the  computations 
end  describee  the  gain  in  fatlgu*  performance  obtained  by  the  reduction 
in  the  gross  allowable  design  screes. 

At  this  time  weight  estimates  were  conducted  for  ell  fatigue- 
critical  primary  structure  on  the  reference  airplane,  end  the  increase 
in  structural  weight  for  e  given  reduction  in  groee  allowable  stress 
was  determined.  With  this  information  sod  the  strsse  increment  versus 
life  increment  date  given  in  Pig.  29,  it  was  possible  to  study  the 
effects  of  the  reliability  analysis  on  design  weight  and  then,  by  raf 
erence  to  Fig.  X),  to  translate  this  into  the  effects  on  payload  or 
range. 

Figure  31  gives  the  relationship  b# tween  reliability  and  weight 
Increase  for  a  structural  component  designed  by  fatigue  cone  id erst  lone, 
whereas  Fig.  32  shows  e  plot  of  the  percent  increase  in  structural  weight 
versus  the  level  of  reliability  of  tbs  weakest  airplanes  in  a  fleet  of 
600  aircraft.  The  initial  weight  of  the  airplane  reflects  currant 
design  practice,  which  includes  the  scatter  factor  of  4  for  fatlgum- 
crltlcal  structure.  The  increase  in  weight  results  directly  from  the 
increase  in  the  scatter  factors  mm  computed  from  Table  XX.  These 
values  were  normalised  to  the  current  factor  of  4,  end  the  resultant 
ratios  represent  the  Increased  life  factors  required  by  the  proposed 
reliability  analysis.  Any  of  the  reference  structural  locations  that 
were  examined  during  this  study  end  that  did  not  mast  these  new  criteria 
had  their  gross  stress  lowered  sufficiently  to  meet  them,  and  the 
Increase  in  structural  weight  wee  subsequently  computed.  It  should  be 
realised  that  aajor  structural  airframe  components  may  have  quite  dif¬ 
ferent  fatigue  performance  capability,  as  some  components  may  be  almost 
entirely  designed  by  fatigue  considerations,  whereas  other  locations 
have  to  seat  static  requirements  primarily,  but  also  have  built-in 
fatigue  capacity.  As  e  coneequenca  of  chess  design  requirements  ell 
the  reference  components  that  ware  analysed  were  not  found  to  be  equally 
penalizing,  end  the  plotted  curves  showing  weight  increase  with  speci¬ 
fied  reliability  (Fig.  32)  reflect  this  aspect  of  aircraft  structural 
design. 

One  other  very  important  assumption  made  for  this  weight  increase 
study  needs  further  clarification.  It  was  considered  essential  that  in 
order  to  test  the  feasibility  of  the  reliability  analysis  an  upper  bound 
condition  was  necessary  to  determine  the  mar  Isa—  penalty  that  could  ba 
incurred  bv  this  proposed  approach.  Consequently,  it  was  asstned  that 


24 


cha  only  way  of  Meeting  the  new  d«al(n  |Ml  waa  by  lowering  chs  gros* 
design  stress  for  all  structural  coaponMti  Chat  war*  both  fatigue 
oriented  and  related  to  the  rtftrnci  local  Iona  chat  had  boon  effect  ad 
by  Ch*  now  analyala.  This  1*  obviously  •  pass  inis  tic  as  sump  t  Ion ,  aa  no 
allowance  was  aada  for  ch*  poaalblllty  of  reaching  the  revised  goal  by 
a  redesign  or  aodlflcation  of  the  existing  detail  structure.  Neverthe¬ 
less,  Fig.  32  dsBonstrates  that,  even  for  this  upper  bound  condition, 
a  weight  penalty  of  3.6X  of  the  primary  structural  weight  was  sufflclant 
to  provlds  s  aadlsn  reliability  on  the  weakest  of  a  flsst  of  600  air¬ 
planes,  and  as  auch  as  932  reliability  was  obtsined  by  s  7X  Increase  in 
the  weight  of  the  prlaary  • true  cure. 

The  lower  curve  on  rig.  32  represents  the  maxim*  weight  penalty 
that  would  be  Incurred  if  the  designer  were  to  specify  reliability  on 
the  second  failure  in  the  fleet.  This  plotted  raault  wa*  based  on  the 
value#  presented  In  Table  X,  which  tabulates  scatter  factors  for  the 
sucond  failure  of  a  fleet.  It  Is  1—s dlttsly  obvious  free  s  comparison 
of  the  weight  plot*  that  if  a  criterion  baaed  on  the  second  failure 
were  acceptable,  a  auch  lower  weight  penalty  would  be  Incurred,  for 
example,  it  can  be  saan  that  the  weight  increment  chat  would  provide 
median  reliability  of  the  weakest  neaber  in  600  airplanes  weald  also 
result  In  an  80X  level  of  reliability  of  the  second  weakest.  Putting 
this  another  way,  *  design  erltarlon  baaed  on  the  second  failure  of  the 
fleet  would  require  about  two-thirds  of  the  weight  Increment  necessary 
for  the  seas  level  of  reliability  of  the  first  failure  In  the  fleet. 

The  plots  in  Fig.  30  give  che  percent  decrease  in  payload  and  ranga 
respective] y  versus  percent  Increase  in  structural  weight.  These  plots 
were  used  to  develop  Figs.  33  and  34,  which  show  the  decrease  in  payload 
or  range  versus  degree  of  reliability  of  the  weakest  airplanes  in  e  fleet 
of  600.  As  In  Fig.  32  the  upper  end  lover  curve*  represent  the  penal cy 
Incurred  by  specif lection  of  the  first  or  second  failures  respectively 
end,  as  before,  an  effective  lowering  of  the  required  additional  weight 
is  demonstrated.  The  maximum  payload  condition  used  In  this  analysis 
amounts  to  approximately  27. SX  of  the  gross  weight  of  the  reference  air¬ 
plane.  It  is  obvious  that  the  plotted  results  of  payload  versus  relia¬ 
bility  are  Intrinsically  dependent  on  this  figure  and  that  other  air- 
planes  of  greater  or  lower  payload  capability  would  be  affected  either 
leas  or  considerably  more  than  the  reference  airplane.  However,  although 
this  argument  la  undoubtedly  true,  it  should  be  rnmamfrarad  that  the  par¬ 
ticular  airplane  was  chosen  because  of  the  availability  of  specific  data 
and  wlehout  prior  knowledge  of  the  outcome. 

It  must  be  remiphaaiaed  that  chess  studies  are  limited  to  slue  ini* 
structures  only,  as  the  reliability  analysis  can  account  for  only  this 
material  at  che  present  time.  However,  the  large  majority  of  the  air¬ 
plane  structure  is  composed  of  this  material  and  consequently  it  is 
believed  that  the  results  presented  are  representative  oc  the  entire 
airplane  and  serve  to  demonstrate  the  feasibility  of  the  proposed  relia¬ 
bility  plan. 
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To  further  Illustrate  the  application  of  tha  rallablllty  plan,  It 
‘a  assumed  that : 

a  Alrplano  average  Ufa  *  40,000  hours 
than,  andar  praaant  cr It aria, 
aafa  Ufa  -  40,000/4  -  10,000  hours 
a  Scatter  factor  from  rallablllty  estimate  obtained  f ran  Tab  la  IX: 

S  •  1.32  (for  ona  taat  l 
°f  datall)  I 

SH  -  4.96  (for  N  -  600)  |  S  -  7.2 

Sg  -  1.102  (mean  Ufa)  J 

a  Haan  Ufa  of  weakest  •  40,000/7.2  ■  1,600  houra 
a  Ratio  of  Ufa  -  10,000/1,600  -  1.79 

That  la,  the  Weibull-distrlbut ion-baaed  rallablllty  criterion 
requires  a  79X  increase  In  life  to  compare  with  tha  currently  specified 
value  of  10,000  houra.  From  Pig.  29,  a  79X  Ufa  inprovsnanc  requires 
a  15X  reduction  in  stress,  which  translates  into  a  3.6X  Increase  In 
structural  walght,  as  shown  In  Pig.  32. 

in  other  words  a  3.6X  increase  in  structural  weight  would  be  auf- 
ficlant  to  provide  enough  reliability  that  the  aean  performance  of  the 
weakest  of  a  fleet  of  600  airplanes  would  equal  or  exceed  10,000  hours. 
Furthermore,  this  lncreaae  in  structural  walght  translates  into  a  2.8X 
decrease  In  payload  or  a  3. OX  decrease  In  range  for  tha  type  of  airplane 
under  consideration.  Moat  Importantly,  it  ehould  be  reenpheelxed  that  the 
nuabar  of  exposed  detail*  (l.e.  fleet  else)  aa  wall  aa  tha  number  of 
test  details  demonstrating  tha  fatigue  performance  level  are  tha  sig¬ 
nificant  e laments  of  this  reliability  approach. 

6.  RELIABILITY  ANALYSIS  BASED  ON  THE  LOG-NORMAL  DISTRIBUTION 


a.  Results  of  Data  Analyses 

Historically  tha  log-normal  distribution  hat  been  acctptod  as  a 
reasonable  distribution  function  for  fatigue  life  even  at  small  prob¬ 
abilities  of  failure.  Consequently,  it  is  a  very  familiar  tool  within 
the  aircraft  industry  and  aa  such  has  bsan  ths  basis  of  many  past  and 
present  assessment  techniques.  Furthermore,  a  great  deal  is  known  about 
estimating  and  establishing  confidence  bounds  for  the  two  normal  param¬ 
eters,  especially  for  uncensored  samples.  For  these  reasons  it  was 
decided  that  the  log-normal  model  merited  investigation  as  a  possible 
basis  for  a  reliability  analysis. 

Ths  approach  used  parallels  the  work  described  in  the  preceding 
discussion  on  the  Welbull  distribution,  except  that  tha  tvo-ordered- 
fal lure-estimation  procedure  was  not  extended  to  the  log-normal  model. 
This  does  not  imply  that  the  work  accomplished  with  this  distribution 
model  was  inferior  to  that  performed  with  the  Welbull  modal  but  simply 
that  the  timely  application  of  order  statistics  was  not  possible  for 
estimating  the  shape  parameter  of  the  parent  population.  As  a  rasult 
this  was  obtained  by  the  HLE  procedure  (Ref.  26)  after  censoring  of 
dats  that  contained  high-time  outliers. 


26 


The  rMulti  of  che  point  MtlMtM  of  the  shape  and  scale  parameters 
sr*  given  In  Appendix  II,  end  the  eat  imatea  of  the  weighted  imi|e 
va .use  of  ehepe  pereeeter  ere  eueeerlsed  In  Table  XII.  Theee  tabulated 
results  deeonetrate  trenda  similar  to  those  obeerved  for  tha  Welbull 
nodal,  namely,  an  lnproveeent  In  che  magnitude  of  the  answers  with  the 
censoring  of  the  high-time  out  Here.  Th  la  effect  la  illustrated  In 
Fi«a.  33  through  3d.  Theee  are  plota  of  the  percentage  frequency  and 
cumulative  frequency  distributions  ot  the  unbiased  shape  parameter  a 
for  the  samples  of  else  two  to  five  specimens  that  do  not  contain  high- 
tine  outliers.  Cumulative  frequency  plota  for  data  that  have  not  been 
censored  for  hlgh-tiae  outliers  era  con pared  on  these  figures  and  can 
be  observed  to  exhibit  greeter  scatter. 

Figures  39  through  42  are  plota  of  the  anpirlcal  and  theoretical 
s  amplify  distributions  for  the  nane  groups  of  alas  two  to  five  specimens . 
These  thaoratical  curves  have  been  generated  from  the  weighted  average 
value  of  all  tha  qualified  data,  that  la,  cr  -  0.133.  (dee  Table  XII.) 

An  interval  estimate  of  this  point  value  was  calculated  fr c-*  assumptions 
of  normality  and  was  found  to  be  ?5  -  0.14.  The  normal  aa sumption 

was  justiflsd,  for  It  was  known  that  ths  astlmate  followed  a  X*~ 
distribution,  which  reduced  to  normality  because  of  tha  very  large  data 
sample.  It  la  worth  mentioning  that  as  a  result  of  the  weighting  tech¬ 
niques  used  In  this  arslyaia,  tha  1,121  groups  of  data  used  in  estimating 
tha  above  values  of  tha  shape  parameter  era  theoretically  the  equivalent 
of  a  teat  smapla  containing  in  excaes  of  3,300  data  pointa.  Curves  rep¬ 
resenting  tha  93X  confidence  upper  bound  value  have  been  Included  In 
Figs.  39  throiq|h  42.  These  figures  damonetrate  that  the  anpirlcal  and 
theoretical  values  are  similarly  distributed  but  that  tha  theoretical 
values  are  displaced  toward  slightly  higher  scatter  except  at  the  upper 
percentiles . 

Figurea  39  through  42  show  that  th*  sample  svaragee  wars  higher  than 
the  interval  estimate  given  by  all  the  qualified  date,  even  though  only 
13X  of  the  data  was  distributed  to  the  right  of  the  values  predicted  by 
theory.  However,  these  data  result  in  increasing  the  estimated  weighted 
average  value  of  the  shape  parameter,  especially  as  tha  satiMt ion  pro¬ 
cedure  used  for  tha  log -no real  model  coa pules  the  variance  rathar  than 
the  shape  parameter.  It  la  believed  that  the  direct  estimation  of  tha 
unbiased  shape  parameter  rathar  than  the  traditional  procedure  of  taking 
che  square  root  of  the  unbiased  variance  would  result  In  a  alight  lower¬ 
ing  in  the  value  of  the  average  shape  parameter  and  consequently  an 
improvement  In  the  proximity  of  the  empirical  end  theoretical  distribu¬ 
tion  plots.  This  expected  trend  can  be  demonstrated  by  referring  to 
Figs.  43  and  44,  which  are  the  plota  of  the  Welbull  shape  parameters 
obtained  by  the  ML£  procedure.  A  comparison  of  theee  graphs  with 
Figs.  40  and  41  demonstrates  similar  bahavior,  axcept  that  tha  Welbull 
empirical  and  theoretical  curves  are  in  closer  proximity  then  those  for 
the  log -normal  case.  This  explanation  does  not  preclude  the  additional 
possibility  that  part  or  all  of  the  improved  fit  is  inherent  to  tha 
Welbull  model  itself. 


In  keeping  with  th«  procedures  described  earlier  for  ths  Wethull 

(odd,  scatter  factors  were  computed  based  on  the  Interval  esc  lasts 
^  •  0,1*.  As  before,  factors  were  generated  for  the  Influence  of 

sample  siae,  degree  of  reliability,  and  fleet  else  factor,  and  are  tab¬ 
ulated  In  Tables  XIII  and  XIV  to  illustrate  the  relative  Importance  of 
the  three  variables.  These  factors  are  noted  to  ba  of  relatively 
equivalent  weight,  unlike  the  Velbull  nodal  In  which  only  a  snail 
penalty  can  be  attributed  to  the  lack  of  an  Infinitely  large  test  aeaiple. 
It  la  also  obvious  that  the  factors  are  generally  of  a  lower  order  than 
the  equivalent  Welbull  values  and  that  they  do  not  Increase  as  rapidly 
with  expanded  fleet  sixes  or  higher  reliability. 


b  Example  Problems 

The  examples  worked  out  on  pages  21, 22, and  23  say  now  be  reevaluated 
using  aeauMptlona  of  log-normality. 


The  first  exasple  concerned  e  trelllng-edge  discontinuity  detail. 
From  full-acele  structural  taet:  Time  to  first  crack  -  18,1>00  hours 
and  low-tlna  failure  In  268  inspected  airplanes  -  2,720  hours.  Froa 
Table  XIII: 

Test  sample  slxe  -  1  failure 
Reliability  of  aadian  life  of  waakaat 
Flaet  alas  of  268  alrplanaa 


Sn  -1.70 
1.0 
2.50 


than  aadian  Ufa  of  waakeat  of  fleat 
and  currant  as t hod  computes  safe  Ufa 


S*  - 

-  18,500/4.25  -  4,350  hours 

-  18,500/4.0  -  4,600  hours 


The  second  exeaple  concerned  e  critical  skln/doubler  detail.  Froa 
full-acala  structural  taet,  tine  to  crack  -  7,500  hours  and  lcw-tlna 
failure  In  500  airplanes  •  655  hours.  Froa  Table  XIII: 

Teat  sample  slxe  •  1  failure  Sn  -  1.70 

Reliability  of  nedien  life  of  weakest  Sr  -  1.00 

Fleet  else  of  500  airplanes  Sns  -  2.65 

then  asdlan  Ilfs  of  weakest  of  fleat  -  7,500/1.70  x  1.0  x  2.65 

-  1,670  hours 

and  current  method  computes  safe  life  -  7,500/4.0  »  1,875  hours 

It  la  again  apparent  that  despite  e  change  in  the  assumed  distribu¬ 
tion  to  the  Log-normal  model,  the  results  still  demonstrate  the  lower 
order  of  reliability  that  can  be  expected  by  application  of  the  arbitrary 
scatter  factor  of  4,  that  is,  lees  than  50X  reliability  on  the  weakest 
members  of  s  large  fleet. 


One  additional  item  naeds  further  discussion,  namely,  the  technique 
used  in  establishing  the  seen  life  of  the  respective  details.  A  compari¬ 
son  with  the  aerlier  examples  on  peg  as  21, 22,  and  23  shows  that  in  the  former 
case  account  wee  taken  of  the  felled  end  uniailed  starboard  and  port  taet 
wings,  whereas  in  the  log-normal  case  the  wing  was  considered  as  one  total 
unit.  This  was  a  necessity  at  this  tlae,  for  although  tools  are  available 
for  obtaining  point  estimates  from  censored  data  of  limited  sample  slxe 
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their  h 1 «m  and  variance  ara  not  known  and  hanca  It  1*  not  currently  poe- 
albLe  to  determine  a  suitable  interval  aatiaata.  However ,  the  teat  vaa 
considered  aa  a  eonpiate  aaatpla  and  by  adopting  the  aaaa  criterion  for 
the  fleet-  that  la,  each  airplane  haa  one  coaplata  wing — it  waa  poaaible 
to  work  the  problaai. 


29 


SECTION  IV 


MATHEMATICAL  RESULTS  FOR  SHAPE-PARAMETER  ESTIMATION 

The  mathematical  results  necessary  to  obtain  a  numerical  estimate 
of  aluminum  alloy  structural  fatigue  "scatter”  are  listed  and  briefly 
described  in  this  section.  With  reference  to  the  schematic  (Fig.  2) 
of  the  reliability  analysis,  this  section  corresponds  to  the  second 
stage  of  the  plan.  Unlike  the  third,  fourth,  and  fifth  stages,  this 
part  of  the  analysis,  once  completed,  does  not  have  to  be  repeated  for 
subsequent  applications. 

1.  STATEMENT  OF  THE  PROBLEM 

From  a  large  number  of  small  samples  of  applicable  fatigue  data 
obtain  a  precise  one-sided  confidence  interval  estimate  of  the  distri¬ 
bution  shape  parameter  for  both  the  Weibull  ard  log-normal  time-to- 
f allure  models. 

2.  GENERAL  SOLUTION 

It  must  first  be  assumed  that  the  k  selected  samples  of  data  come 
from  populations  with  different  scale  parameters  but  equal  shape  param¬ 
eters  when  applying  either  the  Weibull  or  log-normal  model.  The  error 
in  this  assumption  may  be  evaluated  at  the  conclusion  of  the  exercise  by 
comparing  the  theoretical  sampling  distributions  of  the  shape-parameter 
estimator  with  the  observed  distributions  of  the  data  estimates. 

Because  of  the  above  assumption  the  estimator?  of  the  shape  param¬ 
eter  used  must  be  invariant  estimators,  that  is,  independent  of  the 
scale-parameter  values.  Since  the  large  number  of  estimates  of  the  vari¬ 
ous  data  groups  or  samples  can  be  averaged  to  obtain  one  comparatively 
accurate  estimate,  it  is  obviously  advisable  to  work  with  unbiased 
estimators . 


In  addition  to  these  reasons,  invariant  unbiased  shape  parameter 
estimators  will  be  used  so  that  the  following  important  theorem, 2  stated 
without  proof,  can  be  applied  to  provide  interval  estimates. 

Theorem:  Let  be  an  unbiased  invariant  estimate  of  the  true  param¬ 
eter  x  from  the  ith  of  k  data  groups.  If  the  variance  of 
*i  is  equal  to  X2,  then  the  estimate 


*■ .  jp.  >1Aii 

lx  1/Q‘ 

is  an  unbiased  estimate  of  X  with  variance 


(IV-1) 


2 


Refer  to  Mann  (Ref.  17, 


p.  64.1)  for  a  different  iisage  of  this  theorem. 
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(IV-2) 


V*r  ($')  -  - - 

S  ^ 

i-i  1 

Equation  CIV-1)  is  being  used  to  obtain  point  estimates  of  the 
shape  parameters  from  k  selected  subsets.  By  combining  Eqs.  (IV-1) 
and  (IV-2)  and  using  the  Tchebycheff  inequality  (since  the  distribu¬ 
tion  of  V  is  not  generally  known)  we  obtain  precise  bounds  for  X 
given  by 


P  (- 


±_ 
1  +  C 
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*  X  * 


for  C 
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i  1  where 
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d/Qj) 


1 

2 


(IV-3) 


(IV-4) 


Equation  (IV-3)  will  be  used  in  conjunction  with  all  shape-param¬ 
eter  "subestimates"  to  pinpoint  an  accurate,  but  still  conservative, 
value  of  the  shape  parameters  for  use  in  the  reliability  assessment 
plan. 

3.  SHAPE-PARAMETER  ESTIMATORS  TO  EVALUATE  STRUCTURAL  FATIGUE  SCATTER 
Three  methods  for  estimating  fatigue  scatter  have  been  employed  In 
conjunction  with  the  Weibull  model.  These  are: 

a)  The  two-order-statistic  estimator  in  which  only  the  two  lowest 
failure  times  in  each  qualified  group  of  data  are  used  to 
estimate  <. 

b)  The  maximum- likelihood  estimator  with  censoring  in  which  all 
failure  times  In  qualified  complete  samples  are  used  except  those 
high  times  judged,  by  the  previously  specified  procedures,  to  be 
outliers 

c)  The  maximum- likelihood  estimator  without  censoring  in  which  all 
failure  times  in  complete  samples  are  used  to  estimate  k. 

Estimators  of  the  type  corresponding  to  the  categories  in  b  and  c 
were  aiso  employed  to  estimate  a,  the  log -normal  shape  parameter.  It 

probably  would  have  been  enlightening  •  _v  ..c  and  a  vsaful  two- 

erdtt -statist ic  estimator  fcr  the  log-normal  model  if  time  had 
permitted . 


It  is  emphasized  that  in  all  cases  the  bias,  variance,  and  sam¬ 
pling  distributions  of  these  invariant  estimators  were  required  and 
vs  re  calculated  by  means  of  existing  theory,  when  possible,  or  by  mean s 
of  Moote-Carlo  simulation  adopted  for  this  study. 

4.  WEIBULL  TWO-ORDERED-STATISTIC  ESTIMATOR 

All  statistical  analyses  in  this  project  have  been  directed  toward 
estimating  the  earliest  failure  times  in  a  fleet  and,  accordingly,  it  is 
pertinent  to  consider  estimating  the  scatter  itself  from  only  the 
earliest  two  failures  in  each  of  the  selected  groups  of  data.  Since 
it  has  been  demonstrated  (Sec.  Ill)  that  the  Veibull  estimates  are 
most  affected  by  relatively  unimportant  high-time  outliers,  it  is 
especially  relevant  to  estimate  the  Veibull  shape  a  from  the  first  and 
second  time- to-f allure  order  statistics.  This  approach  should  be  nearly 
Insensitive  to  the  effects  of  high-time  outliers  that  cannot  be  explained 
by  the  Veibull  model.  The  final  weighted  two-failure  estimate  will  be 
compared  with  the  weighted  likelihood  estimate  using  all  the  data  both 
with  and  without  the  arbitrary  procedure  for  censoring  outliers  des¬ 
cribed  in  Sec.  III. 

As  explained  before,  a  best  unbiased  invariant  estimator  of  <  is 
sought  so  that  Eqa.  (IV-1)  and  (IV-3)  may  be  used  to  obtain  a  weighted 
average.  Following  is  an  outline  of  the  derivations,  developed  by 
Mann  (16) ,  of : 

a)  The  estimator  itself  [Eq.  (IV-12)] 

b)  The  sampling  distribution  of  the  estimator  [Eq.  (IV-13)] 

c)  The  variance  of  this  sampling  distribution  [Eq.  (IV-14)] 

Let  X  and  X_  be  the  first-  and  second-ordered  observations 
l,n  z,n 

respectively  from  a  sample  of  n  independent  Weibull  random  variables, 
each  with  frequency  function 


where  a,  8 


Given 

statistics 


(n-2)  ! 


fl- 1  r  Q  -i 

fOc)  m  j  (j)  exp  l-(^)  J  for  X  >  0 
>  0  are  the  parameters  of  the  distribution. 


Eq.  (IV-5) ,  the  joint  density  of  the  first- 
can  be  derived  and  is  given  by 


(IV-5) 


and  second-order 


(IV- 6) 
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where:  0  <  X,  *  x  <  « 

l,n  2  ,n 

The  density  of  the  scale  invariant  acatlstlc  Z  ■  X,  /X,  can  be 
derived  f  roe  Eq .  (IV-6)  as  n  z*n  *»n 


,  n(n-l)  a  i 

*(z)  ’  r~~n-m5n  1 

It  fellows  that  the  distribution  of  Z  la 

n 


for  a  s»  1 


(IV- 7) 


CU>  -  P  (7, 


a 

-  n  f  — iSTilJ-  *  ! - — 

*'l  [l  +'  ( n- 1 )  u]  Z  1  +  (n~l 


(IV-8) 


•»Ince  has  scale  Invariance,  it  follows  that  any  estimator  of 
<  -  1/a  that  is  a  function  of  Zn  alone  is  also  invariant.  Hann  (16) 
his  derived  the  maxlamg- likelihood  estimator  of  <  fro*  the  first  two 
of  n  failures  and  found  It  to  be  expressible  in  the  form 


A 

k  ■>  u  in  Z 
n  n 


(IV-9) 


where  the  coefficient  uR  depends  only  on  n.  This  led  Hann  to  consider 
all  estimators  of  the  same  general  form  as  in  Eq.  (IV-9),  namely ,  the 
product  of  a  constant  depending  on  n  and  tn  Zn.  Among  these  estimators 
there  is  only  one  unbiased  invariant  estimator  «c ,  and  it  may  be 
expressed  as 


^  •  0  tn  Z 
n  n 


(IV- 10) 


so  that  the  coefficient  bn  is  giver,  by 


n  E  (tn  Z  ) 

n 


It  should  be  noted  now  that  there  is  no  difference  between  this 
unbiased  two-ordered  statistic  estimator  ic  and  the  unbiased  likelihood 
estimator  V,  which  will  be  introduced  later,  provided  that  k  la  calcu¬ 
lated  at  the  second  observed  failure  time.  However,  there  will 
generally  be  a  small  numerical  difference  since  the  to-be-discussed, 
coefficients  arc  determined  approximately  by  computer  simulation, 
whereas  the  coefficients  bn  can  be  calculated  exactly  from  the  results 
that  follow. 


Rl  (*)--£  &=*±L  for  x  el 

J”1  j" 

Equations  (IV-l«s!)  and  (IV-I4)  have  been  used  with  the  weighted 
average-type  estimators  in  Eqa.  (IV-1.)  and  (IV-3)  to  determine  « . 
Equation  (IV-1 3) ,  the  theoretical  distribution  of  the  subest vaatoi  , 
has  been  compared  with  the  observed  distribution  of  the  data  group 
estimates  in  figs.  15  through  18.  The  resulting  close  fit  between 
theory  and  observation  attests  to  small  error  in  the  Initial  assumption 
regarding  the  existence  of  a  unique  value  of  «c  that  controls  the  scatter 
of  fatigue  lives. 
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5.  WEIBULL  MAXIMUM-LIKELIHOOD  ESTIMATOR 

Prior  to  application  of  the  tvo-ordered-etatlatic  type  of  tmtinator 
to  mluiod.se  the  problem  of  high-fat lgue-11 fe  outlier*,  an  estimation 
procedure  was  sought  that  would  best  account  for  all  fatigue  life 
observations  in  each  sample.  The  KLE  was  a  suitable  choice  since,  for 
general  application,  it  provides  either  a  beat  estimate,  or  nearly  so, 
and  for  specific  application  to  Weibull  shape  estimation  it  outperforms 
all  investigated  alternatives.  (Refer  to  Ref.  18.) 

On  this  basis,  when  high-life  outliers  do  not  occur,  which  has 
apparently  been  the  case  in  tests  on  large  mult idetai led  components, 
the  KLE  is  recommended  for  pinpointing  the  shape-parameter  value. 
However,  since  the  vast  majority  of  the  data  consists  of  small,  mono¬ 
lithic  specimens  of  the  type  that  contains  a  minority  of  definite  high 
outliers,  the  cvo-ordered-statlstic  estimator  was  relied  on  in  this  work 
to  obtain  the  estimate  of  <  most  relevant  to  the  lower  portion  of  the 
fatigue  life  distribution. 

As  discussed  before,  this  two-ordered-statistic  estimator  is  simply 
a  special  case  of  the  MLE  where  all  fatigue  life  observations  in  a  sam¬ 
ple  have  been  ignored  except  for  the  two  lowest.  In  view  of  this 
relationship  and  because  the  problem  of  high-tiae  outliers  nay  be 
coped  with  in  auy  number  of  other  ways  (such  as  by  looking  at  all 
observations  except  those  judged  to  be  outliers)  ,  it  was  considered 
worthwhile  to  discuss  some  MLE  results  obtained  during  this  task  and 
prior  to  application  of  the  order- stat 1st ic  approach. 

a.  Weibull  MLE  Equations 

The  derivation  of  the  estimators  given  below  for  a,  6  -*  0  was 
first  given  by  Cohen  (13)  and,  for  completeness,  is  Included  In 
Appendix  I. 
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i-J 
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£  *n  Xt 

i-1 


(IV-15) 


Equation  (IV-15)  la  not  explicit  in  x  and  aunt  be  solved  by  an 
iterative  procedure  for  each  nev  data  sampling.  The  correct  value  of 
6  aay  than  be  substituted  into  Eq,  (IV-16)  to  obtain  the  MLE  of  the 
Weibull  scale,  censed  6.  This  estimate  0  is  to  be  used  only  whan  x  is 
assumed  to  be  unknown.  For  the  case  to  be  considered  in  the  next  sec¬ 
tion,  that  of  known  r,  the  following  equation  is  applicable: 


A 

8 


(IV- 17) 


A  V  A 

b.  Sampling  Distribution  of  Weibull  MLE  x.  8,  and  6 

It  Is  believed  that,  except  for  theAcase  of  ti(  *  2  studied  by 
Mann  (16),  the  distributions  of  the  MLE  x  and  $  of  Weibull  distribution 
parameters  have  never  been  calculated.  However,  the  distribution  of  6 
la  known  precisely  for  samples  chat  are  either  complete  or  censored 
by  failure  time,  and  its  formulation  is  included  in  Appendix  I. 

The  key  to  the  concist  description  of  all  three  sampling  distri¬ 
butions  is  to  express  the  estimator  within  an  invariant  or  parameter- 
free  staciatic  and  to  find  the  distribution  of  that  statistic. 

For  the  previously  solved  case  of  8  it  is  shown  in  Appendix  I 
that  the  known  distribution  of  the  statistic 

«  1/k 

W-  <£)  ( IV-18) 

D 


is  Independent  of  true  values  of  parameters  a  and  8.  Likewise,  it  is 
shown  that  the  unknown  marginal  distributions  of 

A 

U  -  -  -  T  ( IV-19 ) 

a  x 


and  of 


«  l/* 
V 


t.IV-20) 


do  not  depend  on  a  and  8.  A  method  is  then  outlined  in  Appendix  I  for 
finding  these  distributions  approximately  through  Honte-Carlo  computer 
simulation.  The  results  of  the  work,  where  2,0cJ  Independent  random 
samples  were  generated  per  distribution,  are  presented  in  Tables  III 
and  IV  and  in  Figs.  4  and  5.  Notice  the  excellent  agreement  between 
the  empirical  and  the  exactly  known  distributions  of  U  for  samnles  of 
two  observation,!. 
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c.  Bias  and  Variance  of  the  ^atimator  k  ■  1/a 

The  bias  and  variance  of  tc ,  in  terns  of  *  ,  have  been  evaluated 
for  several  complete  saeple  sites  free  the  empirical  distributions 
of  1/U,  The  results  are  given  below  in  a  form  easily  combined  with 
Eqs.  (iV-1)  and  (1V-3)  to  obtain  an  overall  estimate  of  k  from  all 
qualified  complete  samples; 


Complete 


Bias 


Variance 


Sample  Size  n 

Factor  B 

n 

Factor  < 

2 

1.73 

0.71 

3 

1.37 

0.35 

a 

1.25 

0.22 

5 

1.187 

0.164 

10 

1.088 

0.073 

20 

1.047 

0.033 

where  E(B  r)  -  k  and  Var  (B  k)  -  Q  r2. 

n  n  n 

The  workable  suj^estiuate  is  then  formed  by  multiplying  the 
likelihood  estimate  k  by  the  coefficient  to  obtain 

hi  : 

n 

an  invariant  unbiased  estimate  of  k. 


6.  LOG-NORMAL  MAXIMUM-LIKELIHOOD  ESTIMATORS 

Assume  that  the  logarithm  XL  of  the  life  span  is  a  random  variable 
with  a  normal  distribution,  where  u,  o  >  0  are  the  two  parameters  of 
the  distribution.  Cohen  (26)  has  derived  the  maxisium-likelihood  esti¬ 
mators  u  and  a  of  u  and  a  for  incomplete  samples,  as  given  by 
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(IV-21) 


where  is  the  estimated  hazard  function  or  failure  rate  at  time  Z^ 
of  one  of  the  ng  test  terminations  for  reports  other  than  failure. 
Thus , 
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♦  U^ 

-  i~Tc  jT 


♦  (ei)  «  /  1  ♦  (t) 


dt 


where : 


G.  A 

i  *  U 

ci  — 


♦  (t)  -  — - exp 

v:2f 


[A] 


end  where  and  s  are  the  seen  and  variance  of  the  log  lives  ot 
specimens  that  tailed;  that  is 


“f  \ 

\  *  E  ‘M  * 

1*1  f 
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£  (X.  -  X,)'  /n 

i-1  Li  1 


Neither  C  nor  a  is  given  explicitly  by  Eq.  (IV-21)  for  n_  *  1,  so 
a  oust  be  solved  oy  an  iterative  procedure.  For  complete  samples 
(n^  ■  0),  Eq.  (IV-21)  reduces  to 

U  -  XL  ( IV— 2 2 ) 

and 

a2  -  s2  (IV-23) 


The  bias,  variance,  and  sampling  distributions  of  the  estimator 
in  Eq.  (IV-23)  are  well  known  (for  example,  refer  to  Ref.  27)  and  are 
listed  below; 


a. 


Bias  of 


a2 

a 


so  that 


E(o2) 


n-1  2 

n 


2 

a 


in  which 


( IV-24 ) 
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tharef ora, 

2 

Var  <»2> '  t£t  (a2)  (IV-25) 

For  th«  case  of  nail,  incomplete  samples,  ch«  bias  and  variance 
of  che  estimators  in  Fq.  (IV-21)  are  not  known,  but  additional  computer 
simulation  studies  could  be  attempted  to  provide  these  values.  At  this 
tine  only  the  complete  samples  and  the  censored  samples  with  more  than 
50  observed  failure  times  are  used  to  estimate  a. 


The  resultant  0.95  confidence  interval  value  o  >0.14  was  used  In 
all  log-nonaal  applications  that  follow  and,  somewhat  arbitrarily,  for 
application  of  the  Tchebvcheff  inequality  to  gage  the  effect  of  assuming 
no  failure  model. 
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SECTION  V 


INTERVAL  ESTIMATION  OF  CERTIFIABLE  LIFE  FOR  THE  FLEET 


All  the  estimation  and  probability  thaory  corresponding  to  the  last 
thrae  etapa  of  the  analysis  is  discussed  in  this  section.  Referring 
again  to  Fig.  2,  these  steps  are: 

Step  3:  Obtaining  a  stoat  likely  estimate  of  Che  distribution  of 

fatigue  lives. 

Step  4:  Moving  b»ck  Co  a  high  confidence  bound  estimate  of  the 

fatigue  life  distribution. 

Step  5:  Mapping  the  distribution  of  step  4  into  an  order  statistical 

distribution  from  which  a  life  may  be  certified  for  the 
safety  of  all  or  moBt  of  the  exposed  fleet. 

To  facilitate  these  steps,  which  must  be  repeated  with  each  appli¬ 
cation  of  the  plan,  all  relevant  equations  have  been  translated  into 
simple  relations  froa  which  scatter-factor  values  may  be  easily  calcu¬ 
lated.  The  total  scatter  factors  are  henceforth  defined  asthe  ratio 
of  one  of  the  best  estimators  of  central  life  (either  §,  1CP ,  or  10*) 
and  the  calculated  safe  life 

1.  STATEMENT  OF  THE  PROBLEM 

From  one  saaq>le  of  structural  fatigue  data,  use  a  selected  t lot¬ 
to-failure  model  with  known  shape  parameter  to  specify  a  valid  proce¬ 
dure  for  estimating,  to  any  desired  confidence  level,  either  (1)  the 
distribution  of  the  Mc^-ordered  failure  time  in  &  fleet  of  size  N  or 
(2)  simply  the  distribution  of  an  arbitrary  (unordered)  failure  time. 
Choose  one  of  these  distributions  at  a  desired  fr&ctile  to  certify  a 
design  life. 

2.  GENERAL  METHOD  OF  SOLUTION 

To  obtain  as  much  accuracy  as  possible  or,  more  precisely,  to 
obtain  uniformly  minimum  (mean-squared)  error,  maximum-likelihood  point 
estimators  of  the  distribution  scale  parameters  will  be  used.  The 
sampling  distributions  of  these  (unbiased)  estimators  are,  for  the 
most  part,  known  so  that  interval  estimates  of  the  scale  parameters 
may  be  calculated  exactly  fros  the  point  estimates. 

The  known  shape  parameter  and  the  scale-parameter  interval  esti¬ 
mate  completely  determine  die  confidence  interval  bound  on  the  distri¬ 
bution  of  uaordered  times  to  failure.  Since  the  shape  parameter  is 
assumed  to  be  known,  this  bound  has  the  same  failure  distribution 
(Ueibull  or  log-normal)  as  the  point  estimate,  except  that  it  is  scaled 
backward  in  time  by  a  constant  multiplying  factor. 

The  fractiles  of  the  order-statistic's  distribution  may  then  be 
calculated  precisely  by  means  of  the  well-known  cumulative  binomial 
formula.  This  completes  the  theoretical  application.  A  fleet  life  can 
then  be  -estimated  that  corresponds  to  administrative  specifications. 
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3.  BEST  ESTIMATE  OF  FATIGUE-LIFE  DISTRIBUTION 

With  che  shape-par ammtar  value  ash used  known,  we  depend  on  the 
fuil-ecale  fatigue  teat  to  obtain  the  most  likely  value  of  the  parent 
population's  central  tendency  or  scale  parameter.  For  this  reason  it 
is  recognised  that  the  occurrence  of  high-time  outliers  in  this  testing 
could  lead  to  overly  high  safe-life  estimates.  Fortunately,  as  has 
been  observed  in  the  data,  large,  complex,  multicomponent  specimens 
are  not  characterized  by  outliers  to  anywhere  near  the  extent  that 
simple  monolithic  specimens  are. 

This  will  complete  Che  Job  of  obtaining  the  fatigue  life  distri¬ 
bution  most  representative  of  the  data  for  the  given  two-parameter 
failure  model. 

a.  Point  Estimate  of  Weibull  Characteristic  Life 

As  given  in  Sec.  IV  [Eq.  (IV— 17)  ],  the  MLE  of  0,  with  a  ■  1/k 
known,  is 


This  estimate  is  workable,  in  the  sense  that  its  Important  prop¬ 
erties  arc  known,  only  for  those  small  samples  that  are  complete  o 
censored  by  failure  time.  Considering  the  case  of  single-stage  ce,  >r- 
ing  at  the  nfth  failure  time  only,  Eq.  (V-l)  becomes 


t 


(n-nf) 


1/a 


(V-2) 


Equation  (V-2)  gives  an  unbiased'  uniformly  minimum  variance 
estimator  for  which  the  sampling  distribution  is  known. 


b.  Point  Estimate  of  Log-Normal  Median  Life 

For  small  samples  (n^  < SO) ,  only  the  case  of  uncensored  samples 
will  be  considered  here.  The  relevant  MLE  is  simply  the  data  log- 
average 


A 

y 


y 

i-i  108  Y: 


( V-  3) 


« 


This  scale  parameter  estimator  is  free  of  bias  and  possesses  minimum 
variance  of  expected  loss  (mean-squared  error). 


^Actually  E  (t)  *  0  and  E  (0a)  ■  3° 
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c.  Point  Intimate  of  Single  I’arawtor  :  Ti  hebvi  hot  f  Inequality 

For  this  case  of  no  assumed  failure  model,  wc<  will  choose  to 
estimate  the  average  of  logarithmic  failure  times  rt  ■  log  Y.  Fills 
Is  consistent  with  the  arbitrary  use  of  o  ■  0.14  as  the  known  value 
of  the  standard  deviation  of  logarithms.  The  arbitrariness  arises 
from  the  fact  that  this  particular  value  was  influenced  by  previous 
assumptions  of  log  normality  and  here  we  are  applying  it  as  a  "known" 
quantity  in  a  distribution-free  limit  theorem. 

4.  LOWER  BOUND  INTERVAL  ESTIMATE  OF  FATIGUE-LI FK  DISTRIBUTION 

It  is  in  this  step  that  knowledge  of  the  scatter-controlling  shape 
parameter  yields  the  greatest  dividends.  For  one  thing,  the  analysis 
is  greatly  simplified;  but  far  more  importantly,  the  lack  of  a  sampling 
of  many  full-scale  structural  fatigue  life  observations  becomes  a  small 
rather  than  an  insurmountable  problem.  The  sampling  distributions  of 
the  Weibull  §  and  log-normal  0  estimators  are  given  below  in  terms  of 
a  probability  statement  of  confidence  level  and  in  terms  of  a  scatter 
factor  s  *  accounting  for  finite  sample  size  alone.  The  same  task 

is  done  with  the  data  leg-average  estimator  n,  using  distribution-free 
methods . 

a.  Exact  Lower  Confidence  Bounds  for  Weibull  scale  Parameter 

Appendix  I  contains  the  proof  of  a  theorem  Important  for  our 
applications;  namely,  that  the  statistic  2nf  W  •»  2nf  (8/B)a  has  the 
chl-squared  distribution  with  2rif  degrees  of  freedom.  It  is  there¬ 
fore  possible  to  express  an  interval  estimate  of  6  as 


*;  <2nf) 


(V-4) 


where  (2n^)  is  the  y-fractile  of  the  chi-squared  variate  with  2n^ 
degrees  of  freedom  and  where  B^  satisfies 


P(8y  <  B)  -  y 


The  scatter  factor  S  can  then  be  defined  that  shows  the  penalty 

n 

paid  to  gain  confldenca  y  from  a  finite  sample  size.  This  factor  is 
expressed  as 


(V-5) 
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b.  Exact  Lower  Conft  'nee  Round*  for  Log-Normal  Smlt*  Parameter 

One  of  the  moat  widely  uaed  statistics  te  ^ ,  which,  for  the 

case  of  complete  samples,  has  standard  normal  distribution.  It  Is 
therefore  possible  to  obtain  an  exact  lower  hound  of  p  as 


C  ♦  k 


(V-b) 


where  the  standard  normal  deviate  k  is  a  negative  number  when  y  Is 
greater  than  0.5.  Y 

Equation  (V-6)  leads  directly  to  the  statement  of  confidence 


p  (uy  *  W>  -  Y 


(V-7) 


since  fro.*  properties  of  the  normal  distribution  we  can  say 

P  (p  »  P  ♦  k  0/  JH  )  •  y 

y  v 


(V— 8) 


Since  p  and  p  are  in  terms  of  log  life,  we  define  S  as 
Y  — *  n 


S  -  10 

n 


(p-p  )  ..  -k 

Y  *  10  yjn 


(V-9) 


c.  Precise  Confidence  Bounds,  No  Distribution  Assumed 

If  absolutely  nothing  is  known  or  assumed  about  the  distribution 
of  fatigue  lives  except  the  .alue  of  its  (logarithmic)  variance,  one 
(weak)  tool  that  My  be  applied  to  bound  n  ts  the  distribution-! ree 
Tchebycheff  inequality.  If  the  variance  of  log  lives  is  (j2f  the 
variance  of  n  is  known  to  be  o^/n,  regardless  of  the  distribution.  We 
may  then  write,  for  any  number,  6*1, 


P(  In  -  n|  s 


)  »  1 


(v-10) 


This  allows  the  construction  of  a  lower  bound,  n^,  given  by 


*  n  - 


(V-ll) 


/  n  (  1-y j 


which  satisfies 


P(n^  t  n)  *  y 


(V-12) 
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W«  cun  evaluate  the  scatter  (actor  l.y  means  of  Eq .  ( V  -  A 1 )  as 


(V-t 3) 

5.  FACTORS  FOR  LOWER  BOUND  ESTIMATE  OF  A  CERTIFIABLE  LIFE 

At  this  stage,  the  final  part  of  the  analysis,  all  estimation  work 
has  been  concluded.  In  all  that  follows,  the  lower  bound  Interval 
estimate  Py(y)  is  treated  mathematically  as  If  It  were  the  true 
fatigue-life  distribution  7(y).  Of  course,  the  resultant  safe-life 
calculation  la  only  an  aatlaata  at  the  same  chosen  level  of  confidence  y. 

It  is  assumed  here  that  the  following  specification,  recommended 
in  Sec.  II,  is  made  by  customer  and/or  vendor  management  for  a  particu¬ 
lar  component. 

With  probability  t,  the  safe  life  that  la  being  eatiaiated 

at  confidence  level  y  shall  be  exceeded  by  at  least  N-M+l  of 

the  N  exposed  structural  components. 

The  first  cask  is  to  transform  F  -  1-R,  the  failure  probability 

of  the  weakest  of  N  Identical  components,  into  an  equivalent  value  of 
F,  the  failure  probability  of  any  randomly  selected  component.  This 
may  be  done  precisely  through  the  cumulative  binomial  formula,  whose 
derivation  appears  in  Ref.  28. 

The  probability  of  exactly  M  failures  and  the  probability  of  at 
least  H  failures  (equal  to  the  failure  probability  of  the  "M™  weakest" 
member)  are  given,  respactlvely,  by 

P(T  -  M)  -  (JJ)  F*  (1-F)  N_M  (V— 14) 

P(T*M)  -  £  (J)  F1  (1-F)  N~1  (V— 15) 

i-M 

A  fora  of  Eq.  (V-15)  that  is  easy  to  calculate  for  saall  values 
of  M  and  large  N  is 

P(T*M)  F  -  1-  £  (?)  F1  (1-F)""1  (V-16) 

M.N  i-0  i 

This  aquation,  as  nantionad  before,  was  plotted  for  several  values 
of  M  and  M  in  Fig.  3. 

Strictly  speaking,  Eq.  (V-16)  is  not  tbs  cumulative  binomial  dis¬ 
tribution  as  implied  but  alaply  one  minus  same.  It  can  always  be  used 
to  calculate  the  "general"  reliability,  associated  with  a  random  fleet 
member ,  that  corresponds  to  the  specified  reliability  for  the  Mth 
failure  among  N  members  for  all  integer  values;  thus. 


10 


(n 


V  -  ia(s/n  (1  -  y)  ^ 
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1  <  M  «  N  •  • 

However,  since  Che  design  parameter  advocated  Is  usually  associated 
with  the  first-  or  second-order  statistic,  attention  may  be  restricted 
to  the  special  cases: 

(1)  M  -  l,  N»1 

where  Eq.  (V-16)  reduces  to 

F  ■  Fj  -  1— <1— F)N  (V-17) 

or,  in  terms  of  reliabilities,  to 

R  -  RN  (V-18) 

(2)  rf  -  2,  N*2 

where  Eq.  (V-16)  reduces  to 

F  -  F2;K  -  1-(1-F)N  -  HFd-F)*-1  (V-19) 

or,  in  terms  of  reliabilities,  to 

R  -  NRN_1  1 1-R  (-^p)  |  fV— 20) 


a.  Estimate  of  Safe  Life:  Weibull  Model 

As  shown  previously,  the  lower  bound  interval  estimate  of  reli¬ 
ability  R  as  a  function  of  component  age  y  is  given  by 

R  (y)  -  1-F  (y)  -  exp  |-(y/H)°|  (V-21) 

Assume  now  that  the  minimum  acceptable  value  of  R_haa  been  calcu¬ 
lated,  as  a  function  of  specified  values  of  M,  N,  and  R,  by  means  of 
Eq.  (V-16).  Then  we  need  only  solve  for  y  in  Eq.  (V-21)  to  obtain  the 
"inverse  reliability  function"  and,  therefore,  the  safe  life 


T  •  Y(R)  -  8  [in  (1/R)]l/* 


(V-22) 


The  scatter  factor  accounting  for  both  the  fleet  exposure  and 
the  specified  fleet  reliability  R  is  simply 


s„.i  -|t.  (1/«)| 


-1/a 


(V-23) 
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For  the  case  where  the  entire  fleet  (weakest  member)  is  being 
protected,  we  have  M  ■  1  and 


R 


R  1/N 


So  that  safe  life  and  scatter  factor  are  given,  respectively,  by 


and 


-  In  (1/R) 

N 


—  in  (1/R) 

N 


1/n 


-1/a 


( V-  2  4 ) 


( V-  2  5 ) 


For  thia  case,  the  safe  life  and  the  total  scatter  factor  S  can 
be  written  as 


and 


V 


tn(l/R)  11/a 

*>f)  J 


(V-26) 


S  -  S 


2  / n  1/a 

«,  »vi 

tn  (1/R) J 


( V— 2  7) 


b.  Estimate  of  Safe  Life:  Log-Normal  Modal 

The  development  for  the  case  of  the  reliability  being  specified 
in  terms  of  the  time  to  first  failure  is  no  more  difficult  for  this 
model  than  for  the  Weibull  model.  The  resultant  safe  life  and  scatter 
factors  are  listed  below: 


?„  -  iot“  +  vl  .  10  ''*)”] 


1/N 

where  R  ■  R 


(V-28) 


(V-29) 


(V-30) 
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c.  Kstimate  of  Safe  Life:  So  Modal 

We  again  make  use  of  the  generally  applicable  Tchebychet f  limit, 
theorem  to  establish  the  reliability  function  for  this  ease  as 

2 

R(y)*l  -  — *——2 - -  t  log  y  <n  -  o  (V-31) 

(- •  log  y) 


and  for  known  a. 

The  safe-life  and  scatter-factor  equations  can  then  be  derived  as 


An  - 


yr  ■ 10  jr,i/u 


) 


1-R 


or 


V  '  10 


( V—  32 ) 


(V—  33) 


Jn  ~(l-y)  Jl  -~R 


1/N 


and 


SR  " 


10 


10 


[jrnr17"] 


( V-  34 ) 


and,  therefore,  the  total  scatter  factor  is 


S  *  Sn  '  SR  * 


10' 


10 


[>/n  (1-y)  \1  -  R1/NJ 


(V-35) 


To  facilitate  numerical  comparisons  and  safe-life  computations, 
scatter  factors  corresponding  to  various  values  of  all  relevant  param¬ 
eters  were  computed  and  presented  In  Tables  IX,  X,  XIII, XIV,  and  XV 
and  in  Fig.  45  through  50.  In  the  tabulations  it  was  arbitrarily 
decided  to  subdivide  S_  into  two  factors,  S-  and  S  ,  such  that  for 
fixed  H,  *  R  ^ 


S 


R 


SR 


This  was  done  to  enable  a  comparison  of  the  factor  Sg  required 
for  reliability  in  just  a  one-member  fleet  and  of  the  factor  S 

K 

required  to  account  for  the  size  N  of  the  exposed  fleet.  The_  factor 
&NSf  as  can  be  inferred  frou  the  double  subscripts,  depends  on  R  as 
well  as  on  N.  This  condition  arises  because  the  order  statistic's 
distribution  is  usually  of  a  differant  fora  than  the  parent  distribution. 
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so  that  one  staple  (actor  (scalar  t rar.s lat  ion )  cannot  be  used  to  com¬ 
pute  all  fractllea  of  the  order  statistic's  distribution.  One  exception 
to  this  is  the  Weibull  model  for  M  ■  1,  where  the  distribution  of  times 
to  first  failure  is  also  Weibull  with  the  same  shape  as  the  parent.  In 
this  case  we  take  the  liberty  of  changing  the  notation  Sfjg  to  to 
ayabollre  the  lack  of  dependence  of  upon  R. 

The  parametric  study  of  these  scatter  factors  brought  out  several 
expected  trends.  The  log-normal  model,  with  low  failure  rates  or 
mortality  at  old  age,  predicts  more  chance  for  a  long-lived  test 
specimen  and,  therefore,  higher  values  of  Sn..  The  Weibull  model, 
with  monotonically  and  rapidly  increasing  failure  rate,  leads  to 
scatter  factors  SNg  and  Sg  higher  than  those  of  the  log-normal  model, 
especially  when  the  fleet  is  large  or  the  specified  reliability  is  high. 
Therefore,  for  the  typically  large  (N>100)  fleets  of  aircraft  In  opera¬ 
tion,  the  Weibull  model  would  lead  :o  a  higher  total  scatter  factor 
S  and  a  lower  safe-life  estimate  chan  ths  log-normal  would  for  a  given 
apecif 1 cat  ion. 
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Srs  I! ON  '.I 

CONCLUSIONS  AND  RKCOMMENI'ATIONS 


1.  CONCLUSIONS 

A  study  of  the  available  options  for  formulating  a  reliability 
plan  based  on  the  collected  aluminum  fatigue  test  data  has  induced 
the  following  conclusions: 

a)  The  most  feasible  method  of  using  the  Weibull,  or  log-norrrval ,  titne- 
to-failure  model  as  a  basis  for  aircraft  reliability  assessment 

is  to  first,  calculate  and  fix  shape-parameter  values  from  a  large 
mass  of  applicable  fatigue  data  and  then  to  base  all  subsequent 
analytical  work  on  these  "known"  values. 

b)  Considering  that  the  applicable  fatigue  performance  data  com¬ 
prised  samples  from: 

•  all  the  investigated  aluminum  alloys,  which  contained  both  the 
major  2000-  and  7f)00-series  alloys; 

•  all  the  notched  configurations,  including  simple  monolithic 
notched,  the  less  simple  structural  simulators  such  as 
laboratory-type  lap  joii.CE,  butt  joints,  etc.,  and  the  full- 
scale  structural  components  and  structures; 

•  axially  loaded  and  flexurally  loaded  test  programs; 

•  constant-amplitude  and  variable-amplitude  testing  procedures; 

and  ^  5 

•  a  fatigue-life  range  between  10  and  10  cycles, 

it  was  demonstrated  that  the  estimated  shape  parameter  was  not 
influenced  by: 

1)  material  alloy  type 

2)  specimen  geometry 

3)  fatigue  test  macnine 

A)  teat  loading  procedure 

5)  cyclic  life 

c)  Considering  the  very  general  conditions  listed  in  conclusion  b,  it 
is  believed  that  for  most  aluminum  structural  applications  a 
reasonable  value  of  the  Weibull  shape  parameter  is  a  *4.0  and  for 
the  log-normal  case  a  •  0.14. 

d)  Fatigue  scatter  is  greater  for  the  following: 

1)  rotating  beam  testa 

2)  hand  forgings 

3)  compression  testing 

4)  adhesively  bonded  joints 

5)  lov-amplitude/high-life  testing 

e)  Fatigue  scatter  is  lower  for  the  following: 

1)  unnotched  specimens 

2)  ultracarefully  manufactured  test  specimens 

f)  Currently  the  two-ordered-statistic  estimator  7  is  the  preferred 
estimator  for  the  task  of  determining  the  Weibull  shape  parameter 
a  from  fatigue  data,  because: 

1)  <  la  simple  to  calculate,  as  no  iterative-type  solutions  are 

necessary; 
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2)  the  Idea  of  emphasizing  the  lovesr.  failure  times  in  each  data 
group  to  calculate  fatigue  scatter  is  consistent  with  the 
concept  of  specifying  safe  life  as  a  fractile  of  the  distribu¬ 
tion  of  the  lowest  failure  time  in  the  fleet; 

3)  an  estimation  procedure  that  considers  only  the  lowest  two 
failure  times  must  be  nearly  insensitive  to  the  occurrence  of 
hightlae  outliers  that  are  not  predicted  by  the  failure  time 
model;  and 

4)  it  was  demonstrated  that  a  fairly  close  agreement  exists 
between  the  observed  and  theoretical  sapling  distributions 
of  fc,  which  lends  credence  to  the  belief  that  the  fatigue 
performance  of  the  weaker  mosbers  of  a  sample  or  a  fleet  may 
be  adequately  described  by  the  Weibull  model. 

g)  In  the  majority  of  cases,  the  use  of  the  Weibull  distribution 
model  will  lead  to  a  lower  estimate  of  safe  life  than  that 
determined  from  the  log-non  1  model . 

h)  The  Weibull  distribution  is  ^accerized  by  a  failure  rate  that 
increases  monotonical ly  and  rapidly;  consequently,  the  calculated 
Weibull  scatter  factors  are  largely  dependent  on  specified  (high) 
reliability  and  fleet  exposure. 

i)  The  log-normal  distribution,  by  reason  of  a  failure  rate  that 
eventually  decreases  to  zero,  predicts  more  chance  for  possible 
high- time  test  failures  than  does  the  Weibull  model.  The  calcu¬ 
lated  log-normal  scatter  factors  reflect  this  and  demonstrate 
greater  dependence  on  fatigue  test  sample  size  than  the  corre¬ 
sponding  Weibull  values. 

J)  The  results  from  a  reliability  analysis  based  on  either  of  the 

considered  distribution  models  relative  to  the  current  approach  to 
safe-life  estimation  (i.e.  the  arbitrary  scatter  factor  of  4)  shows 
that  the  latter  approach  may  result  in  a  very  low  degree  of  reli¬ 
ability  of  the  weakest  members  of  the  fleet,  and  this  is  emphasized 
by  the  large  size  of  some  of  the  current  operational  fleets. 

k)  The  initial  results  of  the  reliability  analysis  demonstrate 

some  promise  of  correlation  between  service  and  predicted  fatigue 
performance. 

l)  The  application  of  the  proposed  analysis  could  result  in  an 
Improvement  in  fatigue  reliability  at  the  cost  of  some  Increase 
In  structural  weight,  with  subsequent  loss  in  payload  or  range. 

Having  demonstrated  the  applicability  of  the  reliability  pier., 
the  following  recommendations  can  be  made. 

2.  R£C0(MENDAT10NS 

a)  The  analysis  should  be  extended  to  cover  some  other  common  air¬ 
craft  materials ,  such  as  titanium  and  high-strength  steel. 

b)  The  Monte-Carlo  studies  should  be  continued  to  complete  the  listing 
of  bias  and  variance  values  for  shape-parameter  estimates  from 
censored  samples  of  small  size  for  inclusion  into  the  computer 
program  to  provide  the  capability  of  assimilating  these  data  into 
the  overall  weighted  estimate. 

c)  Otnei  distribution  functions  should  be  considered,  especially  those 
of  a  fairly  general  nature  but  not  necessarily  as  Indefinite  as 

the  Tchebycheff  limit  theorem  or  other  distribution-free  assumptions. 
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It  Is  expected  that  the  collected  fatigue  data  .ould  he  used 
as  a  guideline  toward  selection  of  representative  node  1  $ . 

d)  The  possibility  of  devising  and  applying  a  two-ordered-stat ist ic 
estimator  for  the  shape  parameters  of  the  log-normal  and  other 

t ime-to- f ai lure  models  should  be  investigated. 

e)  The  results  of  this  study  should  be  incorporated  into  nreviously 
explored  failure  models  that  combine  the  chance  static  and 
fatigue  modes  of  failure. 
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CURRENT  APPROACH  SCATTER  FACTOR  METHOD 


y 


NOTE  Tho  design  life  is  taken  to  be  v<>  the  inference  seems  to  be  that 
PI Y  >  Y<j)  is  acceptable  although  undefined 


Figure  1  Schematic  Representation  of  Current  Reliability  Plan 


PROBABILITY  OF  FAILURE.  F  1  RELIABILITY 


PROCEDURE 


1  Assume  fatigue  life  is  a  random  variable  with  a  two-parameter  Weibull  distribution 

2  Estimate  the  tfiape  parameter  a  from  ail  applicable  data  and  resard  it  as  fixed  in 
all  that  follows. 

3  Estimate  the  scale  parameter  0  from  full  sale  fatigue  test(s).  thereby  defining  Curve  A 

4  Using  the  small  full  scale  test  sample  only,  obtain  a  lower  interval  estimate  0y  of  0 , 
thereby  defining  Curve  B  at  the  confidence  level  V 

5  Specify  or  compute  a  desired  reliability  R  of  an  arbitrary  fleet  member  Use  Curve  B 
to  obtain  the  corresponding  certifiable  life  Y R 


Specification  of  Desired  Reliability 

The  value  of  R.  the  minimum  allowable  reliability  of  •  randomly  chosen  fleet  member, 
may  be  arbitrarily  spaci  jt  may  be  bend  on  and  calculated  precisely  from  a 
specification  of  the  iy  W  of  the  M**1  weakest  of  N  fleet  members. 


Figure  2  Outline  end  Schemetic  Repreeentetion  of  Proposed  Refiebdity  Plen 


PROBABILITY  OF  AT  LEAST  M  FAILURES  IN  A  FLEET  OF  N.  F 


RELIABILITY  OF  A  RANDOMLY  CHOSEN  FLEET  MEMBER.  R  *  1  F 


1/106  1/10*  1/10.000  1/1,000  1/100  1/10  1.0 
FAILURE  PROBABILITY  OF  A  RANDOMLY  CHOSEN  FLEET  MEMBER.  F 

Figure  3.  Distribution  of  Severe/  Order  Statistics 
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RELIABILITY  OF  THE  M*  WEAKEST  IN  A  FLEET  OF  N.  RM>N  •  1  FM;N 
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Figure  5.  Empirical  Distribution  of  the  MLE  of  the  Waibull  Seal* 
Paramatar  0  for  Complete  Samp/as  of  Sira  n  •  m  •  2.  3, 

4,  5,  10,  and  20  (Waibull  Shape  Assumed  To  Be  Unknown ) 
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CUMULATIVE  FREQUENCY 


Figure  6.  Theoretic*/  Distribution  of  the  ML  £  of  the  Weibu/I  Scste 

Parameter  3  for  Complete  Samples  IffeibuU  Shape  Assumed 
To  Be  Known! 


CUMULATIVE  FREQUENCY 


EXPECTED  LOSS  AMO  VARIANCE  Of  lot  /4> 


Figure  7.  A  Mature  of  Sampling  Error  of  the  MLE  of  the  Weibu/I  Shape  From 
Complete  Semples 


ALUMINUM  FATIGUE  TEST  DATA:  TOTAL  OF  EXAMINED  DATA 


A0N3nO3H3  3AIJLVinwnO  MO  A0N3fTO3Hi  3DV1N3DM3M 


f>9 


All  Qualified  Data 


(Lwnittd  to  nmpt«  notched.  iohiU.  and  ftructuraa,  excluding  low  emplitude.  high  -eyeW  data) 


for  Several  Alloys  Several  Specimen  Types 


ALUMINUM  FATIGUE  TEST  DATA 
(Limited  to  simple  notched.  joints,  and  structures. 

excluding  low amplitude.  high -cycle  data)  (Sample  categories  o»  retected  data) 


Estimates  of  the  Weibull  Stupe  Parameter  Estimates  of  the  Weibull  Shape  Parameter 

tor  Different  Types  of  Loading  for  Various  Unacceptable  Data 


UNBIASED  ESTIMATE  OF 


Figure  14.  Observed  Variation  of  Fatigue  Scatter  With  Cycfic  Life  for  1,174  Constant-Amplitude 
Samples  of  Mixed  Sites 
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ALUMINUM  FATIGUE  TEST  DATA 


i 


A3N3fX)3Md  3AIJLViriwnD  HO  ADN3n03Hd  3DVXN33W3J 


Figure  15.  Comparison  of  the  Theoretics/ and  Observed  Figure  16.  Comparison  of  the  Theoretical  and  Observed 
Distributions  of  Estimates  of  tha  Weibull  Distributions  of  Estimates  of  the  Weibu/I 

Shape  Parameter  for  Only  Qualified  Data  Shape  Parameter ,  Using  First  Two  Ordered 

of  Sample  Sire  *  2  Failures  of  Only  Qualified  Data  of  Sample  Site 
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Figure  l 7  Comparison  of  the  Theoretical  and  Figure  18  Comparison  of  the  Theoretical  and 

Observed  Distributions  of  Estimates  Observed  Distributions  of  Estimates 

of  the  Weibull  Shape  Parameter,  Using  of  the  Wetbull  Shape  Parameter  Using 

First  Two-Ordered  Failures  of  Only  First  Two-Ordered  Failures  of  Only 

Qualified  Data  of  Sample  Sui  =  4  Qualified  Data  of  Sample  Size  -  5 


ALUMINUM  FATIGUE  TEST  OATA  TOTAL  OF  EXAMINEO  DATA 


A3N30O3H3  3A  11*101*00  HO  ADN3003H3  3DV1N3DH9H 


<>r» 


ALUMINUM  FATIGUE  TEST  DATA:  TOTAL  OF  EXAMINEO  DATA 
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ALUMINUM  FATIGUE  TEST  DATA:  TOTAL  OF  EXAMINED  DATA 
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■FLIGHT  LENGTH  (HOURS)  - 


0  T«ai  ml  TaMoH 

(5)  dtmb  «o  Cnmt  Altitude 

(5)  Outo  m  Crum  Altitude 

($)  Dnaoant  (Gear  Down) 

0  LmttngandTaaiin 


Figure  25  Typical  Flight  Profile  for  a  Military  Jat  Tanker  TTrarmport  Airplane 
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Jat  Tankwr /Transport  Airplane 
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NUMSER  OF  CYCLES  TO  FIRST  CRACK 

Figure  28.  Typical  Family  of  S-N  Diagrams  tor  a  Structural  Component  on  a  Military 
Jet  Tanker  /Transport  Airplane 
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D€ CREASE  IN  GROSS  DESIGN  STRESS  IPERCEKT) 


% 


Figure  29.  Rdetionahip  of  Fatigue  Lift  With  Street  for  the  Fatigue-Cnttcei 

Aluminum  Structure  of  the  Reference  ktiiitery  Tanker fT \ rentport  Airplane 
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I  MCA EASE  IM  A I  AFLAME  STRUCTURAL  ^  INCREASE  IN  STRUCTURAL  WEIGHT 

WEIGHT  (PERCENT)  1  OF  FATIGUE -CRITICAL  COMPONENTS  (PERCENT) 


v  31.  Relationship  of  Structure!  Weight  of  Fatigum-C fittest  Component  With  Dmgrm  of 
Reliability  of  Weakest  Components 


Figure  32.  Reletionthip  of  Primary  Structural  Weight  With  Dagne  of  Reliability  of 
Weakest  Members  in  Fleet.  Military  Tanker  /Transport  Airplane 
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aluminum  fatigue  test  data 
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Figure  «9  Influence  of  Fleet  Sin  on  Scan*  F acton  Required  ro  Provide  SO %  RoHobMry 
for  the  First  end  Second  Fleet  Feduree,  Weibuil  Model 
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Figure  50.  Influence  of  Fleet  Sin  on  Scatter  Factor*  Required  To  Fmeide  50 %  Reliability 
for  tha  First  end  Second  F hat  Failures,  Log-Normal  Modal 
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Simulated  Examples  To  Illustrate  the  Effects  of  Isolated  Long-Life  Specimens 


Table  in.  Seven  Empirical  Distributions  of  the  Parameter- Free  Statistic 
1/U  =  a/ a 


NOTE:  Tabulated  values  are  (l/U).  where  P  U/U  > L.f)  (1/U)|  1  -  > 
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•For  each  sample  size,  l/U  waa  calculated  for  each  of  2,  000  independent  group*  of 
generated  exponential  (o  -1)  variates.  The  empirical  distributions  were  taken 
to  be  the  ordered  arrays  of  these  estimates. 

b  Vs  lues  in  parentheses  are  obtained  from  theory. 
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Table  IV.  Seven  Empirical  Marginal  Distribution*  of  the  Parameter- Free 
Statistic  V  =  (%/(J  ]r 


NOTE:  Tabulated  values  are  v,  where  P  (V<v)  »  V 


S.  Sample 

siz.en 

Coefficient''"''^^ 
of  confidence 

2 

3 

4 

5 

10 

20 

Consort  <1 
3  of  5 

0.  01 

io"1** 

0.  0005 

0.  03 

0.  10 

0.  39 

0. 34 

iu‘5 

0.  10 

0.01 

0.  19 

0.31 

0.  44 

0.  02 

0.  72 

<».  05 

0.  25 

0.  23 

0.  43 

0.  3S 

0.  04 

0.  79 

0.  94 

0.  25 

o. 

0.  01 

0.  97 

0.  97 

0.  93 

I.  00 

0  99 

0.  09 

0.  75 

2.  OS 

1.  04 

1  55 

1.41 

1.20 

1.  17 

1. 24 

0.  5)0 

11.  0 

3.0 

2.  70 

2.  17 

1.00 

1.39 

1 .  SO 

0.  03 

70.  0 

0.  4 

4.  09 

3.  1 

1.  S3 

1.53 

2.  0 

0.  O') 

3.  000. 0 

17.  0 

7.  7 

5.  2 

2.  21 

l .  75 

4.2 

0.  Of) 

10‘ 

SO.  0 

17.  0 

9.  0 

2.  7 

i.  99 

9.  .. 

Smallest  of  2.  000* 

.o'30 

io'° 

5x1 0“fi 

0.  0004  0.  137 

0. 39 

10'10 

Mean 

10"“ 

io‘* 

3.  1 

1. 39 

1. 09 

1.03 

1. 30 

Largest  <>f  2.000* 

to2* 

10° 

79.  0 

100.  0 

7  4 

2.  OS 

400  0 

Variance  of 
flistrilwtion 

(lo4*) 

no10) 

<90) 

(10) 

(0.24) 

(0.079X100) 

*For  each  sample  size,  V  was  calculated  for  each  of  2, 000  independent  groups  of 
generated  exponential  {0r=  1)  variates.  The  empirical  distributions  were  taken  to 
be  the  ordered  arrays  of  these  estimates. 
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Table  V  S<‘im  Ibei-retieal  Distnltulions  >>f  the  Pa ra nu* ter  -  Free  Statist m 

U  .  B  *1  >° 


NOTES  1  Tabulated  values  are  w,  where  P (W <  w)  V 
2.  From  theory.  FfW)  =  1 

Var  (W)  ^  l/nf 


Cneffieient 

of 

ton>nicnn  > 


Numt.w  r  of 
failure  "f)servati<>ns 
sample  n 


(1 

ot 

o  nio 

■>  07 

0  14 

0  26 

0.  41 

u 

(12 

0  020 

0  11 

0  19 

0.  31 

0.  46 

n 

of) 

0  03 

0  1* 

0.27 

0  39 

0.  34 

n 

10 

0  1  1 

0.  2  7 

0  37 

0  49 

0.  S2 

0 

20 

0  22 

0.41 

0  31 

0.  62 

0  73 

(1 

30 

0  30 

0  33 

0  04 

0  73 

0.  81 

0 

30 

o  09 

0  M 

0  *9 

0  93 

0.  97 

0 

70 

1  20 

1  22 

1  21 

1  is 

1.  13 

0 

so 

1  61 

1  30 

1  43 

1.  34 

1  25 

0 

00 

2  30 

1  <M 

i  77 

1  60 

1  42 

0 

03 

3  00 

2  37 

2  10 

1  *3 

1. 57 

0 

3  91 

2  *>2 

2  31 

2.11 

1  75 

0 

4  *>i 

3  32 

2  -0 

2  32 

1  97 

Table  V!  Result?  «>f  Analyse?  Determining  the  Typical  Shape  Parameters  for 
Fatigue  Performance  of  Aluminum  Structures 


Number 

Unbiased  estimates' 

Data  description 

of 

ot  reciprocal 

groriv 

shape  parameter  « 

AU  collected  data 

2.  003 

0.  259 

All  unqualified  data 

705 

0  325 

AU  qualified  data: 

Limited  to  simple  notched,  joints, 

and  structures,  excluding  low- 

amplitude,  high-cycle  data 

1,  290 

0  224 

All  qualified  data: 

2024 

551 

0.  213 

7075 

402: 

0  227 

717S 

125 

0  240 

Monolithic  notched 

444 

0.  240 

Laboratory  structures  simulator 

557 

0  200 

Full-scale  structures 

297 

0  243 

Constant  amplitude 

1,  1T4 

0.  229 

Variable  amplitude 

124 

0.  177 

|  Bonded 

h3 

0.  419 

•  !  Low  amplitude,  high  life 

lit 

0.  419 

*  F’nnotched 

179 

o.  no 

3 

TO  to  iO  cycles 

50 

0.  240 

10^  to  10*  cycles 

160 

0.  215 

10*  to  6  x  10*  cycles 

434 

0.  204 

0  x  10*  to  4  x  10^  cycles 

444 

0.  259 

4  x  10^  to  106  cycles 

- 

70 

0.  227 

•Rejected  data 
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intimate  from  first 

two -ordered  failures  MI,F -censored  M LE -uncensored 


Table  VII.  Comparison  of  Results  Between  Uncensored  Fatigue  Data  and  Data 

Censored  for  High- Time  Outliers 


Size  of 
groups 

(  n  specimens'! 

Wei  bull  model 

— 

Number 

of 

groups 

Unbiased  estimate 
of  reciprocal 
shape  parameter  a 

420 

0.  253 

3 

328 

0.  301 

4 

234 

0.  286 

5 

145 

0.  380 

Weighted 

value  of 

1.250 

0.  300 

all  groups 

2 

420 

0  253 

3 

279 

0  211 

4 

193 

0.  201 

5 

117 

0.270 

Weighted 

value  of 

1,  119 

0.242 

all  groups 

2 

420 

0.253 

3 

330 

0.229 

4 

236 

0.200 

5 

151 

0.227 

Weighted 

value  of 

1,298 

0.224 

all  groups 

-  -  - 

Table  VM.  Comparison  of  Reaulta  From  Total  Collected  Fatigue  Data  ami 
QuaUfied  Structural- Equivalent- Type  Data 


Sample 
size  n 

Number 

of 

samples 

tm  i 

« 

Obaerved 
Var  « 

Theoretical  * 
Va  r  7^ 
(for  «  7) 

. 

Theoretical* 
Var  7 

(for  «  0.  2241 

tj  «s 

O 

524 

0  293 

0.  1056 

0  0611 

0. 0357 

rjS  S3 

3 

431 

0  304 

0.  1526 

0.  0756 

0  0410 

jC  ±  S  <- 

f-  -  «  •» 

4 

494 

0.  215 

0  0855 

0.  0401 

0.  0435 

S  "2 

5 

196 

0  273 

0.  1271 

0  0667 

0  0449 

■g  «  3 

a-s 

2 

420 

m 

0.  0585 

0. 0456 

0  0357 

3 

330 

Kssa 

0  0537 

0  0429 

0  0410 

5  3  «  ■ 

«■'“  3 

4 

236 

0  200 

0.  0391 

0.  0347 

0.  0435 

5 

151 

0  227 

0.  0841 

0.  0461 

0. 0449 

•In  theory,  Var?-Cn«2  !See  Eq.  (IV- 14)1 
where:  C’2  =  0  712 
C3  =  0.  818 

C  .  «  0.  867 

4 

C5  *  0.  895 
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Tested  to  Failure  (Welhul'l  Model) 
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design  or  certifiable  life 


Table  X.  Scatter  Factor*  To  Obtain  Heliability  R  in  the  Second  Weakeat  of  a  Fleet  of  Size  N,  when  n,  Specimens 

Are  Tested  to  Failure  AVeibull  Model)  * 


Tabic  XI — Continued 


97 


5.100  5  5  ...  9,600  I  3,500  to  5,000 


10*’ 


rable  XU.  Comparison  oi  Ret  Between  Lncensored  Fatigue  Data 
and  Data  Censored  for  High -Time  Outliers  (Log-Normal 
Distribution) 


Size  of 
groups 

,'n  specimens) 


Weighted 
value  of 
all  groups 


Weighted 
value  of 
all  groups 


Log-normal  model 

— 

Number  of 
groups 

o’(log  X) 

420 

0.  149 

328 

0. 181 

234 

0.159 

145 

0.208 

1 , 250 

0.168 

420 

0. 149 

278 

0.111 

193 

0. 108 

119 

0.147 

1,121 

0.135 

Table  XIII.  Scatter  Factors  To  Obtain  Reliability  R  in  the  Weakest  of  a  Fleet  of  Sue 
N.  When  n  Specimens  Are  Tested  To  Failure  (Log-Normal  Model! 


>P|  .10  u3is.»| 


Table  XIV  S<  alter  Fac  tors  To  Obtain  Reliability  ft  in  the  Soron<l  Weakest  of  a  Fleet 


% 


Table  XV.  Scatter  Factors  To  Obtain  Reliability  ft  in  the  Weakest  of  a  Fleet 
of  SUe  N,  When  n  Specimen*  A  ire  Tested  to  Failure 
(No  Failure  Model  Assumed ) 


NOTES: 


i. 


2 


Data  based  on  the  Tchebycbeff  Limit  Theorem  in  lieu  of  an  assumed 
failure  model.  The  log -standard  deviation  a  is  assumed  to  equal  0  14, 
the  unknown  log  average  n  is  estimated  at  a  confidence  level  greater 
than  95%. 

lO'j 

Scatter  factor  rr — *  S  S„  Sp 

h 


V  -  average  of  log  Uvea  of  data 
=  safe  life 


Tabulated  values  are  of 


APPENDIX  I 


MATHEMATICAL  DERIVATIONS,  THEOREMS,  AHD  PROOFS 
REQUIRED  FOR  APPLICATION  OF  THE  VEIWULl.  MLE 


INTRODUCTION 

This  appendix  discusses  all  nonobvious  mathematical  statements 
introduced  In  the  text  that  pertain  to  maxi mum- likelihood  estimation  of 
Weibuli  population  parameters.  These  Include: 

a)  A  derivation  of  the  Weibuli  MLE  for  all  samples. 

b)  Proof  that  the  following  statistics  are  independent  of  a  and  ft, 

for  all  complete  and  some  censored  samples, 

1)  U  -  o/u  A 

2)  V  -  (fi/8)a 

3)  U  »  (ft/e>a 

c)  A  brief  description  of  the  computational  procedure  used  to  obtain 
the  empirical  marginal  distributions  of  U  and  V. 

d)  Some  remarks  pertaining  to  the  adequacy  of  the  random-number 
generator  used  to  generate  independent  exponential  variates  needed 
to  Initiate  the  above  computational  procedure. 

e)  A  derivation  of  the  theoretical  distribution  of  W. 


DERIVATIONS 

Let  X  be  a  random  variable  with  a  Weibuli  distribution,  the 
reliabiXlty  of  X  then  im 

R(x)  *  exp(-(x/8)°J  for  x  >  0  vAl-l) 

where  o,8  >  0  are  the  two  parameters  of  the  distribution. 

We  assume  that  a  number  of  components,  say  n,  each  with  ihe 
same  Weibuli  distribution  of  life  are  put  on  test.  However,  what  is 
ob served  at  each  UAsi  is  the  time  the  component  fails  or  the  time  Che 
tes^  is  terminated.  Let  Z  be  a  random  time  at  which  the  test  is 
terminated  for  any  reason  other  than  failure  of  the  component.  Then 
what  we  observe  Is  the  event 

(X  -  x)  n  (X  <  Zj  or  [Z  ■  i]  O  jx  *  Z)  IAI-2) 

We  now  state 

L«mc>a  1:  if  is  a  vector  of  nonnegative  random  variables, 

possibly  dependent  on  "(X^ , . . . ,Xq) ,  which  are  themselves  independently 
and  Identically  distributed  noonegative  randort  variables  with  common 


density  function  f  and  distribution  F,  then  the  likelihood  of'  the 
event 


k  n 

r  (X  <  z  )  IX  -  X  |  n  (X 
i-l  lilt  1-k+l 


-  z , 


v  A 1  —  3  > 


where  k  ( «  n)  is  the  (random)  number  of  failures  observed,  is  of 
the  form 

k  n 

c  TT  f  (X  >  TT  [i  -  F ( z  > j 

i-l  j-k+1  ^ 


and  the  constant  C  depends  on  but  net  on  F. 

Proof:  Let  g  be  the  joint  density  of  (Zp.  ..  .Z^)  given  X  ,...,X  . 
Then  the  probability  of  the  event  specified  in  Eq.  (Al-3)  is 


ui,xt 


i-l 


k 

TT  f  u.  > 

i-i 


1 1— F(z  )  ]g(z. 

i-k+i 


z  ix,, 
n  1 


T> 


dz^v«  » « >d zk 

which  upon  sitapli f ication  shows  that 


.  / 


V*i 


i-l 


,k  ’ 


'&U, 


,  x  )d; 
n 


,dz, 


(AI-d) 


and  thus  for  an  observed  failure  time  fx. . x  ),  ;  does  nut  depend 

on  F. |  !  Ik 

To  simplify  the  notation  in  what  follows  we  let  denote 

the  set  of  complete  observations  of  X  and  » •  • .  ,  t  )  denote  the  set 

of  censored  tests  (i.e.the  set  of  observations  7).  Then  from  Che 
lemma  we  see  that  the  log-likelihood,  except  for  some  constant  inde¬ 
pendent  of  a,  o,  is 


i-i  6 


t  t  j 

+  (a~l)in(~)  -  (~)  I 

tt  o 


n 

I 

k+1 


.V 

c 


(Al-5) 
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and  e  then  obtain 


H. 

J  J 


■k. 


in(~) 


I 


1 


in(— ) 

o 


»  V  . 

7  Z  ^ 

k  i  * 


(AI-6) 


(AI-? 


V 

Thus  the  two  ma xi mum- 1 i ke 1 i hood  estimators  a,  si  of  a  and  d  when 
both  parameters  are  unknown  are  defined  by 


and 


i-i  = 


( AI-8) 


(Al-9) 


These  results  have  been  given  previously  by  Cohen  (13). 
Let  us  define 


v . 
l 


for  i*l 


u 


y 

V  -  (-) 

c 


then  Eqs.  (Al-B)  and  (AI-9)  beconse,  respectively — by  using  the  fact  that 


and  simplifying: 


(AI-10) 
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and 


n 

i 

i  1 


in  y, 


n 


1 


u 


J. 

k 


k 

L  in 

i-i 


(AI-11) 


For  i-l,...,k  we  know  that  y^  i.«  an  independent  observation 
fro*  an  exponential  random  variable  with  unit  mean.  Consider  the 
condition: 


(Y 


K+l 1 


••V 


CY 


{*' 


•V 


(AI-12) 


has  a  distribution  independent  of  a  and  6,  where  K  is  the  random 
number  of  uncensored  items. 


We  can  now  state: 

Theorem  1:  A  necessary  and  sufficient  condition  that  U  and  V,  as 
defined  in  Eqs.  (AI-10)  and  (AI-11),  have  a  distribution  independent  of 
a  and  6  is  that  Eq.  (AI-12)  be  true. 

Proof:  It  is  clear  from  Eqs.  (AI-10)  and  (AI-11)  that  V  and  V  are 

parameter  free  if  Yj,...,Y  are,  which  is  true  iff  condition  (AI-12) 
is  j  |  . 

Note,  for  example,  that  condition  (AI-12)  is  true  if  with  probability 

one 


K  »  k  and  Z.  . , 
k+1 


(AI-13) 


t  h 

where  x,  is  the  k1-  -ordered  observation  out  of  n.  This  is  merely 

%*n  th 

the  case  of  sampling  censored  at  the  k  failure  time. 


Let  us  briefly  consider  the  case  when  Eq,  s.AI-13)  is  true.  Under 
this  condition  we  can  obtain  the  distribution  of  T  and  V  by 
Monte-Carlo  simulation  methods: 

Let  n  *  k  be  given  positive  integers.  Let  >'(i)  »•  •  •  »y(k)  be 

the  first  k-ordered  observation  from  n  independent  machine-generated 
exponential  variates  each  with  unit  mean.  Solve  for  u  in  the  equation 
♦(u)  -  0,  where 
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(u) 


I  y“l)  a,  y(l)  ♦  <n-k)/‘k)  in  y 


(k) 


5>U)  +  (n-k)y(k) 


u 


i  c- 


£‘"  ycn 


Now  compute 


-I  V  Vu  +  i"-k)vu  u 

k  V  yu)  u  •  (k) 


Repeat  thij  procedure  m  times  for  independent  samples  of  y’s, 
and  so  obtain  a  sample  of  size  m  of  (U,V>.  The  sample  size  m  can  be 
made  sufficiently  large  to  determine  the  distribution  of  U,V  to  a 
degree  of  accuracy  determined  by  the  machine  procedure  that  generates  the 
pseudorandom  numbers. 

The  random-number  generators  used  in  Che  simulation  studies 
presented  here  are  of  the  type  called  composite  congruential  generators. 
These  "second  generation"  methods  appear  to  be  better;  that  is,  they 
satisfy  more  stringent  statistical  tests  of  randomness  than  those  of 
the  simple  congruential  generators  used  previously.  In  the  particular 
method  adopted  here,  three  generators  are  mixed  for  the  IBM  3f>0,  each  of 
which  will  produce  a  full  period  of  residues  relatively  prime  to  the 
modulus  2^  (consequently,  these  mixed  generators  will  produce  2^® 
distinct  random  numbers  before  repeating).  This  method  is  presented  in 
detail  and  discussed  by  Marsaglia  and  Bray  (29),  To  obtain  our 
exponential  variates  with  unit  mean,  we  merely  take  the  negative  of 
the  natural  logarithm  of  the  uniform  variates  generated  by  the  mixed 
congruential  method. 

We  now  consider  the  sampling  distribution  of 

Theorem  2:  If  -  is  known,  the  maximum- likelihood  estimate  of  c, 
call  it  is 
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has  a  distribution  Independent  of  .»  and  d  if  condition  (Al-12)  Is 
true.  Moreover,  In  this  case  the  distribution  is  the  same  as  that  01 


i  n 

«-  k  E  h 

1  1 

Proof:  The  proof  is  Immediate. 

In  the  special  case  when  Eq.  (Al-13)  is  satisfied,  then 

k 


w  ■  i  £  Yu>  *  <-k)Y 


(k) 


Y  ,  ^  are  the  first  k-ordered  observations  from  independent 


where  the 
exponential  variates  with  unit  mean . 

k 


Thus : 


U 


r  £  (n-j+1) [Y 

j-1 


-Y  1 

(J)  (j-l)J 


which  by  known  results  (e.g.  Ref,  30,  p.18)  is  the  sum  of  k 
Independent  exponential  variates  with  unit  mean.  By  the  well-known 
fact  (e.g.  Ref.  31)  that  two  times  an  exponential  variate  with  writ  mean 
is  a  chi-square  variate  with  tvo  degrees  of  freedom,  and  the  reproducing 
property  of  chi-square  variables  under  convolution  (Ref.  30,  p.  10> , 
there  follows  the  result  that  2kW  has  a  chi-square  distribution  with  2k 
degrees  of  freedom. 


HO 


APPENDIX  1  l 


LIST  OF  SALIENT  FEATURES  AND  UNBIASED  POINT  ESTIMATES 
OF  POPULATION  PARAMETERS  OF  COLLECTED  FATIGUE  DATA 


1 .  TABULATED  RESULTS 

The  salient  features  of  all  -he  collected  fatigue  performance  data, 
including  the  estimated  population  parameters,  are  tabulated  in  this 
appendix.  The  estimates  of  the  scale  parameters  were  included  to 
complete  the  description  of  the  data  only,  and  no  further  use  is  made  of 
Che  generated  values. 

An  11-diglt  description  number  is  provided  to  catalog  the  variables 
of  specimen  thickness,  material,  grain  direction,  type  of  structure, 
type  of  specimen,  finish,  cvpe  of  loading,  and  type  of  f  allure.  A  code 
is  provided  on  the  following  oage  to  facilitate  usage  of  the  tabulated 
data  descriptions. 

Also  listed  on  the  output  are  the  test  sample  size  and  the  number 
of  samples  tested  to  failure.  Two  additional  columns,  labeled  "size 
used  for  MLE,  log-normal  and  Weibull,"  are  given.  These  columns 
demonstrate  the  number  of  specimens  used  for  estimating  the  shape 
parameter  after  the  high-cine  outliers  have  been  censored  from  the 
sample.  It  can  be  seen  from  the  tabulated  results  that  the  high-time 
outlier  afflicts  only  a  minority  ot  the  data  and  that  a  sample  censored 
for  the  Weibull  model  Is  not  necessarily  censored  for  the  log-normal 
model . 

The  initial  cwo  shape  parameters  are  the  unbiased  point  estimates 
for  the  log-normal  and  Weibull  models  using  the  MIE  on  data  that  had 
been  censored  for  high-time  outliers,  and  the  iasr  column  is  the  unbiased 
point  estimate  of  the  Weibull  shape  parameter  given  by  the  two-ordered- 
statistic  estimator. 

Finally,  a  complete  listing  of  the  data  references  is  presented,  and 
corresponds  with  the  REF  column  of  the  computer  printout. 
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APPENDIX  III 


LISTED  VALUES  OF  FATIGUE-LIFE  OBSERVATIONS  FOR  ALL  COLLECTED  SAMPLES 


This  appendix  rstabulatss  ths  less  and  description  number*  given  in 
Appendix  II  to  provide  a  cross-referencing  capability  between  the  tvo 
tabulation*.  In  addition!  this  tabulation  itemises  all  ths  individual 
observations  from  the  collected  data  and  categorizes  them  as  either  a 
failed  or  a  suspended  item. 
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*1  A*»T»»CY 

The  application  of  reliability  analysis  methods  to  the  estiaatlon  of  probab!* 
aircraft  structural  fatigue  performance  was  Investigated.  Us*  was  made  of  order 
statistics  to  sstablish  ths  scans  of  assessing  the  fatigue  performance  reliability 
of  a  fleet  or  number  of  fatigue-exposed  details.  A  reliability  analysis  plan  for 
application  to  aluelmas  alloy  structural  fatigue  performance  was  devsloped  and 
compared  with  the  current,  fixed  scatter-factor  procedure  for  determining  the 
safe  life  of  a  structural  detail.  Both  the  two-parameter  Welbull  distribution 
and  the  log-normal  distribution  with  empirically  dtfined  shape  parameter  were 
utilized  to  make  the  reliability  plan  tractabla  as  compared  to  e  distribution- 
f'«e  approach.  ! lex imum- likelihood  estimators.  Including  one  which  considers 
only  the  first  two-ordered  failures,  were  employed  to  examine  the  many  variables 
that  might  Influence  fatigue  scatter,  to  qualify  fatigue  date  that  represented 
aluminum  structural  scatter  end  to  tstabllch  shape- parameter  values  that 
typified  structural  fatigue  scatter.  The  sampling  distributions  of  these 
estimate  s  were  required  to  work  the  problem  and  ware  calculated  by  amen*  of 
existing  theory  or  Monte-Carlo  simulation.  More  than  2,000  groups  of  fatigue 
performance  data  were  collected,  analysed,  and  used  to  dsnonstrate  the  feasi¬ 
bility  of  establishing  a  shape-parameter  value.  Baaed  on  this  estimate,  scatter 
factors  have  been  generated  to  account  for  the  penalty  of  limited  input  infor¬ 
mation,  the  degree  of  desired  reliability,  and  the  else  of  the  exposed  fleet. 

Using  these  factors,  the  possible  effects  of  the  reliability  analysis  on  structural 
weight,  payload,  or  range  wees  explored  for  e  jet-engined  military  tanker/ 
transport- type  airplane. 
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